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ABSTRACT 


-'■Ain  determining  policies  for  the  acquisition  and  management  of 
repairable  spares  for  the  Space  Shuttle,  two  objectives  are  paramount. 
First  is  the  optimization  of  some  measure  of  system  performance  such  as 
the  expected  number  of  shuttles  launched  on  time  per  year.  Second,  since 
the  cost  of  a  spares  mix  can  run  into  the  hundreds  of  millions  of  dollars, 
we  would  like  to  minimize  the  cost  of  achieving  a  certain  performance 
level.  Both  requirements  suggest  a  need  for  mathematical  models  of  the 
supply  system. 

The  high  cost,  low  demand  rate  items  found  on  the  shuttle  are  usually 
controlled  via  an  (s-l,s)  inventory  system.  An  (s-l,s)  policy  involves 
sending  an  item  to  a  repair  depot  immediately  upon  failure.  Using  an 
assumed  (s-l,s)  repair  policy,  this  thesis  will  examine  ways  of  choosing 
a  spares  mix  according  to  three  different  mathematical  models  of  system 
performance  .^Developed  by  Muckstadt  [9],  these  models  are  specifically 
adapted  to  the  ^Derating  characteristics  of  the  Space  Shuttle.  Features 
of  these  models  include  a  nonstationary  Poisson  demand  rate  whose  para¬ 
meter  depends  upon  the  pre-launch  maintenance  schedule  and  a  variable 
weight  of  backorders  over  the  pre-launch  cycle.  Items  are  maintained  at 
predetermined  points  in  time,  and  we  expect  more  demands  on  the  supply 
of  repairable  spares  for  an  item  whenever  it  undergoes  maintenance. 

The  interlaunch  cycles  are  probabilistic  replicas  of  one  another  and  so 
form  a  convenient  time  span  over  which  to  evaluate  system  performance. 
Near  the  end  of  such  cycles,  the  backorder  cost  increases  sharply,  and 
so  the  models  allow  for  changing  weights  of  their  respective  objective 
functions. 


Each  model  generates  spares  mixes  at  various  budget  levels,  and 
the  performance  of  each  mix  is  evaluated  and  compared  with  the  performance 
of  more  elementary  models.  The  models  used  for  comparison  include  the 
one  in  use  by  NASA  when  our  study  was  begun  and  a  Lagrange  multiplier 
technique  based  on  backorders.  We  use  these  models  to  demonstrate  that 
nonstationary  demand  rates  are  important  only  for  long  interlaunch 
cycles  and  for  short  repair  times. 

Another  issue  in  the  minimization  of  delayed  shuttles  is  the  shipment 
policy  used  by  the  serviceable  spares  supply  system.  Since  the  shuttles 
will  some  day  operate  out  of  two  geographically  separated  locations, 
planners  have  the  option  of  building  spares  facilities  at  one  or  both  of 
the  sites.  They  must  also  decide  whether  to  initiate  a  lateral  resupply 
capability  in  order  to  allow  base  to  base  shipments  when  desired.  A 
computer  program  based  on  the  need  and  reluctance  formulas  proposed  by 
Miller  [7]  is  implemented  to  investigate  these  questions.  Without  regard 
to  the  costs  of  the  various  shipment  policies,  the  results  indicate 
that  the  best  performance  is  attained  when  both  bases  have  spares  facili¬ 
ties  and  lateral  resupply  capability.  The  next  best  expected  performance 
comes  from  prepositioning  spares  facilities  at  both  bases  without  the 
lateral  resupply  option.  Least  desirable  is  a  system  with  only  one  fully 
equipped  base  and  lateral  resupply  capability. 
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Introduction 


The  Space  Shuttle,  a  space  transportation  system  introduced  in  1981, 
presents  planners  with  several  unique  features  in  the  area  of  spare  parts 
provisioning  for  repairable  items.  Since  the  system  is  in  some  respects 
very  like  an  airplane,  we  might  initially  concern  ourselves  with  adapting 
the  models  used  by  military  and  civilian  air  operations.  Common  elements 
of  all  these  systems  include  a  number  of  very  expensive  components  with 
low  failure  rates,  a  requirement  that  all  components  must  be  operating 
in  order  to  take  off,  and  a  certain  penalty  for  any  delay.  Thus  many 
such  systems  employ  what  is  known  as  an  (s-l,s)  inventory  system,  whereby 
if  the  stock  level  s  of  any  spare  decreases  by  one  unit,  the  failed  unit 
is  immediately  sent  to  a  depot  for  repairs.  Since  the  repair  times  are 
not  nearly  as  variable  as  the  inter-failure  times,  repair  times  often 
are  fixed  while  failures  are  assumed  to  occur  at  a  Poisson  rate. 

If  the  (s-l,s)  policy  is  pursued  for  shuttle  spares  at  a  given  launch 
site,  there  will  be  two  important  departures  from  the  typical  methods  of 
setting  spares  levels.  First,  the  simple  or  compound  Poisson  demand  rate 
for  failures  will  be  replaced  with  a  nonstationary  Poisson  demand  rate 
reflecting  higher  failure  rates  at  certain  stages  in  the  pre-launch 
countdown.  Since  most  models  employ  demand  rates  over  repair  time  as 
a  measure  of  how  many  spares  are  needed,  we  must  find  some  way  of  incor¬ 
porating  the  time-dependent  demand  patterns  into  our  model.  Secondly, 
the  penalty  for  backorders  will  also  depend  upon  time,  due  to  the  nature 
of  the  pre-launch  sequence.  This  is  especially  true  near  the  launch  date, 
when  delays  may  cause  some  very  costly  preparations  to  be  extended  or 
reaccompl ished. 
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A  second  critical  option  available  to  planners  involves  the  projected 
management  of  spares  assets  when  more  than  one  launch  site  is  operating. 
Since  spares  must  be  shipped  to  a  storage  location  upon  completion  of 
repair,  a  good  shipment  rule  should  account  for  the  proximity  to  launch 
at  each  site  as  well  as  the  spares  stock  and  expected  demand  at  each  site. 
More  important  perhaps  is  the  evaluation  of  the  expected  performance  of 
an  initial  investment  in  two  stocking  locations  versus  only  one,  and  of 
the  option  for  a  lateral  resupply  capability  which  would  allow  one  site 
to  ship  spares  to  another  site  when  necessary. 

The  first  problem,  that  of  setting  spares  levels  for  repairable 
items,  is  solved  by  developing  analytical  algorithms  to  optimize  a  per¬ 
formance  measure  subject  to  budget  constraints.  Each  base  will  be  con¬ 
sidered  separately  although  extensions  are  easily  made,  and  the  depot 
will  be  seen  as  having  infinite  capacity.  These  and  other  assumptions 
will  permit  us  to  focus  on  a  single  inter-launch  cycle  and  to  analyze 
both  nonstationary  demands  over  that  cycle  and  increasing  backorder 
penalties  near  the  launch  date.  We  will  present  details  of  the  models' 
implementations  and  examine  the  patterns  of  spares  mixes  selected  by  each. 
These  patterns  will  be  contrasted  with  two  simpler  algorithms  for  setting 
spares  levels.  From  this  comparison  we  will  show  that  as  the  cycle  length 
shortens,  the  rtonstationarity  in  the  demand  rate  may  be  ignored. 

In  addressing  the  decision  of  how  best  to  design  the  spares  trans¬ 
portation  and  location  system,  we  will  consider  two  sites  operating  at 
different  launch  rates,  and  employ  a  simulation  to  compare  the  performance 
of  the  major  types  of  shipment  disciplines.  We  will  continue  to  assume 
nonstationary  demand  rates,  although  we  will  show  that  this  assumption 
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may  be  relaxed  in  certain  cases.  The  simulation  operates  much  as  a  continu¬ 
ous  review  inventory  policy  would  function  in  real  time.  The  number  of 
spares  in  the  entire  system,  however,  is  fixed,  and  so  the  simulation  is 
not  directly  capable  of  suggesting  both  a  good  spares  mix  and  a  correspond¬ 
ing  spares  distribution  system. 

Some  synthesis  of  the  spares  mix  decisions  and  the  transportation 
and  stocking  location  problem  is  appropriate.  Suggestions  as  to  practical 
applications  of  the  methods  will  be  presented,  as  well  as  some  insight 
into  how  they  might  be  used  interactively  as  planning  tools. 


CHAPTER  1 


The  number  of  spares  required  for  a  given  item  used  by  the  shuttle 
depends  largely  upon  the  failure  pattern  it  experiences  and  upon  its 
importance  to  the  timely  launch  of  the  shuttle.  We  make  the  following 
assumptions  about  the  units  in  question  and  their  operating  environment: 

1.  Each  item  undergoes  maintenance  during  one  or  more  predetermined 
periods  prior  to  launch. 

2.  An  ( s — 1 ,s )  inventory  policy  is  followed. 

3.  The  numbers  of  failures  for  the  different  items  are  independent 
and  have  a  nonstationary  Poisson  distribution  whose  parameter 
declines  as  new  items  are  substituted  for  failed  ones. 

4.  Each  launch  cycle  is  of  fixed  length  and  is  a  probabilistic 
replica  of  all  other  cycles. 

5.  The  repair  facility  has  unlimited  capacity:  each  item  i  has 
fixed  repair  time  T..  which  includes  transportation  time 

to  and  from  the  repair  depot. 

6.  An  unfilled  demand  results  in  a  backorder,  and  backorders  are 
more  critical  late  in  the  cycle.  Backorders  do  not  substantially 
increase  the  length  of  a  launch  cycle,  however. 

7.  Since  resupply  time  between  bases  is  small  compared  to  repair  time, 
a  single  location  is  examined. 

These  assumptions  are  discussed  in  detail  in  this  chapter,  and  a 
general  outline  of  shuttle  pre- launch  operations  is  presented.  A  dis¬ 
cussion  of  models  previously  developed  for  use  in  this  problem  environment 
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is  also  included. 

As  a  first  step  in  the  analysis  of  shuttle  operations,  let  us  con¬ 
sider  the  maintenance  cycle  preceding  the  launch  of  an  individual  shuttle. 

Two  components  of  the  space  shuttle,  the  solid  rocket  boosters  and  the  orb'4-- 
er  are  designed  to  be  reused  many  times.  Beginning  some  time  after  the 
shuttle  returns  from  a  space  flight,  a  pre-launch  cycle  takes  place  during 
which  maintenance  crews  ready  the  shuttle  for  launch  according  to  a 
predetermined  schedule.  Detection  of  a  failed  component  immediately 
leads  to  a  demand  on  the  spare  parts  stock  and  to  the  initiation  of 
repair  on  the  failed  component.  The  number  of  failures  experienced 
during  the  cycle  is  directly  related  to  the  length  of  the  flight  just 
completed,  and  the  components'  failure  rates  are  expressed  in  terms  of 
failures  per  flying  hour.  We  may  assume  that  flights  completed  prior  to 
the  most  recent  flight  do  not  contribute  substantially  to  failures  in  the 
present  cycle,  because  items  are  thoroughly  tested  and  maintained  before 
each  launch.  Although  those  items  which  are  replaced  during  a  cycle  have 
a  failure  rate  somewhat  lower  than  those  which  have  undergone  a  space 
flight,  the  difference  should  be  small  and  is  not  of  great  concern.^ 

The  periods  of  increased  maintenance  activity  are  important  in  that 
they  are  often  accompanied  by  an  increased  number  of  failures,  bringing 

1  The  "new  item"  failure  rate  can  appear  in  another  context.  The  shuttle 
is  composed  of  three  main  units;  the  orbiter,  the  external  tanks,  and 
the  solid  rocket  booster  (SRB).  Because  the  solid  rocket  boosters  are  com¬ 
pletely  overhauled  before  the  launch  cycle  begins,  and  because  a  new  exter¬ 
nal  tank  is  used  for  each  launch,  components  of  these  two  units  have  failure 
rates  independent  of  flying  time. 


about  higher  demands  on  the  stock  of  repairable  spares.  This  is  true  for 
several  reasons.  First,  the  item  may  have  failed  during  flight  or 
earlier  in  the  cycle  but  escaped  detection  (or  remained  inaccessible  to 
the  crews)  until  the  crews  actually  had  contact  with  it.  Second,  many 
of  the  tests  performed  to  insure  that  an  item  is  working  properly  place 
higher  than  usual  stress  upon  the  item  and  may  contribute  to  a  failure. 
Lastly,  there  is  the  possibility  that  the  test  equipment  is  not  working 
properly  and  mistakenly  indicates  that  the  item  is  broken.  This  is 
equivalent,  from  a  spares  standpoint,  to  an  actual  breakdown  of  the  item, 
because  the  item  must  be  sent  to  the  depot  for  tests.  The  shipment 
and  testing  time  may  take  nearly  as  much  time  as  a  normal  repair.  The 
increased  failure  rate  during  these  maintenance  periods  has  been  shown 
to  be  significant,  so  that  it  is  necessary  to  treat  the  number  of 
demands  at  a  given  time  as  a  random  variable  with  a  nonstationary  Poisson 
distribution.  This  is  true  regardless  of  whether  the  item  was  on  the 
shuttle  during  flight  or  was  recently  installed. 

A  random  variable  is  said  to  have  a  nonstationary  Poisson  demand 
rate  if  its  failure  rate  varies  as  a  function  of  time.  Its  distribution 
is  specified  by  the  failure  rate  function  XQ(t).  Through  it  we  determine 
the  failure  rate  over  an  interval  of  length  At  as  follows: 

t+At 

X ( t , t+At )  =  /  X  (s)  ds. 
t 

Since  much  of  our  analysis  is  in  terms  of  expected  demand  during  resupply, 
we  define  A(i,j),  the  approximate  lead  time  (resupply  time)  demand  rate 
for  item  i  evaluated  on  day  j,  as  follows: 
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A(t,t+At) 


A(i.j) 


time 


(a) 


(b) 


Figure  1.1  Failure  Rates  and  Lead  Time  Demand  Rates  as  a  Function  of  Time. 

i 

A(i,j)  *  l  A(t,t+1) 
t=g-T.+l 

where  T..  is  the  resupply  time  for  item  i.  Figure  1.1  illustrates  how 
A(t,t+At)  and  A(i,j)  depend  upon  the  particular  time  interval  being 
examined,  where  m  is  the  day  on  which  maintenance  is  performed,  and  L 
is  the  time  between  launches  for  a  single  shuttle.  We  will  refer  to  the 
sharp  increase  in  demands  on  day  m  as  a  demand  spike. 

>  ,i 

We  will  assume  later  that  failure  detection  does  not  lengthen  the 
maintenance  period  devoted  to  other  units  of  the  same  item  type,  and 
that  the  length  of  the  maintenance  cycle  is  always  fixed  at  L,  and  cannot 
be  increased  by  item  backorders.  These  assumptions  enable  us  to  assert 
that  each  cycle  is  identical  to  every  other  cycle  in  demand  patterns  and 


in  length.  Lastly,  we  assume  that  the  cycles  are  continuously  repeated, 
so  that  the  history  of  demands  prior  to  a  given  cycle  is  identical  in 
probability  to  that  of  every  other  cycle. 

A  failed  item  entering  repair  is  returned  to  the  serviceable  spares 
stock  upon  completion  of  repair.  The  repair  time  T  is  a  random  variable 
and  is  independent  of  the  number  of  items  already  in  repair.  For  our 
purposes  it  is  acceptable  to  assume  that  T  is  fixed,  and  so  items  are 
returned  to  serviceable  spares  stock  T  time  units  after  they  are  removed 
from  a  shuttle.  The  maintenance  crews  will  replace  a  failed  unit  as 
quickly  as  possible,  and  if  there  are  no  spares  on  hand  when  a  demand 
occurs,  a  backorder  results.  There  may  be  some  time  period  d(t)  during 
which  no  penalty  is  incurred  on  such  a  backorder,  but  in  general  d(t) 
is  quite  brief  (or  zero)  due  to  precedence  relationships  in  the  mainten¬ 
ance  schedule. 

While  backorders  on  some  days  may  carry  only  a  small  penalty,  those 
occurring  on  other  days,  especially  just  before  a  launch,  may  be  signifi¬ 
cantly  more  detrimental  to  shuttle  operations.  This  is  because  prepara¬ 
tions  which  are  made  just  prior  to  launch  are  often  more  involved  and 
costly  than  earlier  activities.  These  last-minute  activities  range  from 
temperature  control  of  the  liquid  fuel  to  assembly  of  a  launch  control 
team  consisting  of  many  specialists  and  technicians.  They  are  generally 
not  part  of  the  maintenance  cycle  as  we  have  described  it,  but  are  rather 
more  closely  linked  to  the  operational  activity  of  the  shuttle.  A  delay 
during  this  phase  will  be  very  expensive,  even  though  it  is  not  likely  to 
be  long  enough  to  violate  our  earlier  assumption  of  a  fixed  cycle  length. 
Since  cycles  are  expected  to  be  on  the  order  of  a  week  or  more,  while 
backorders  (as  will  be  seen  later)  last  no  more  than  a  day,  backorders 
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will  probably  not  be  the  cause  of  prolonged  launch  delay.  If  w(t)  is 
the  severity  or  weight  of  a  backorder  at  time  t,  we  can  assume  that  w(t) 
increases  with  time  as  shown  in  Figure  1.2. 


Figure  1.2  Severity  of  a  backorder  (w(t)). 


After  time  L,  we  assume  that  the  shuttle  will  be  ready  for  launch.  The 
actual  launch  takes  place  at  some  time  during  a  launch  window,  which  is 
limited  by  such  factors  as  crew  rest  requirements  for  the  astronauts,  ice 
accumulation  on  fuel  tanks,  weather  conditions,  and  daylight  availability 
at  both  the  launch  and  landing  sites.  If  the  countdown  is  interrupted 
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and  the  launch  window  closes  before  the  shuttle  can  launch,  a  minimum  of 
48  hours  must  elapse  before  the  next  launch  attempt,  resulting  in  another 
setup  and  incurring  a  large  penalty  cost.  This  sort  of  postponement  would 
probably  be  due  to  a  system  failure,  not  a  backordered  spare.  Thus  our 
assumptions  about  w(t)  should  be  reasonable.  In  addition,  since  there  is 
only  a  small  possibility  of  an  item  failure  after  all  pre-launch  mainten¬ 
ance  is  completed,  we  may  assume  that  a  shuttle  will  make  no  more  demands 
on  the  supply  system  after  time  L.  If  a  shuttle  should  experience  a  sig¬ 
nificant  ground  delay  after  time  L,  its  mission  will  be  completed  late, 
but  there  is  sufficient  slack  in  the  schedule  to  allow  that  shuttle's 


next  pre-launch  maintenance  to  begin  on  time.  Thus,  variations  in  the 
schedule  due  to  backorders  or  any  other  problem  do  not  disturb  the  pattern 
of  pre-launch  maintenance  cycles  we  have  described.  All  demands  for  a 
given  shuttle  take  place  during  its  maintenance  cycle,  and  each  cycle 
is  identical  in  probability  to  the  next. 

Current  plans  for  the  shuttle  program  call  for  one  or  more  shuttles 
operating  out  of  two  launch  sites.  A  shuttle  will,  in  general,  return 
to  the  base  from  which  it  was  launched,  and  so  at  each  base  the  shuttles 
will  progress  through  pre-launch  maintenance  in  a  certain  order.  There 
will  be  some  overlap  in  the  launch  cycles,  so  that  the  demand  distribution 
for  a  given  item  on  a  given  day  is  the  sum  of  the  demand  rates  for  items 
on  the  shuttles  being  maintained  on  that  day.  The  number  of  demand  spikes 
experienced  during  any  one  period  L  is  exactly  equal  to  the  number  of 
shuttles  based  at  the  launch  site,  assuming  that  items  are  maintained  once 


per  cycle  on  each  vehicle.  The  demand  distribution  for  each  item  is  still 


nonstationary  Poisson,  but  with  a  parameter  A0(t)  = 


M 

l  A. (t) ,  where  x  .(t) 
i=l  1  1 
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is  the  demand  rate  of  the  it*1  shuttle  on  day  t.  This  superposition  of 
demand  rates  is  illustrated  in  Figure  1.3  for  a  two-shuttle  system  with 
no  time  between  pre- launch  cycles. 


x0(t> 


““ ’  one  shuttle 
-  two  shuttles 


Figure  1.3  Superimposed  failure  rates  for  a  two-shuttle  system. 


Here  the  length  of  the  maintenance  cycle  is  constant  but  the  expected 
number  of  demands  per  cycle  is  twice  as  great.  Note  that  if  we  redefine 
a  launch  cycle  to  be  the  time  between  any  two  launches,  we  experience 
successive  periods  of  length  L/2  whose  expected  number  of  demands 

is  equal  to  the  expected  number  of  demands  for  one  shuttle  over  a  cycle 
of  length  L.  If  the  time  between  launches  is  not  fixed,  we  may  not 
redefine  the  cycle  in  this  way,  and  so  we  find  it  convenient  to  assume 


;«**«u.**'  «v*v- 
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that  the  time  between  launches  is  fixed. 

It  is  useful  to  point  out  some  similarities  between  the  system  des¬ 
cribed  above  and  a  system  in  which  there  is  only  one  shuttle  launched 
at  twice  the  original  rate.  If  the  launch  cycle  is  defined  as  the  time 
between  launches  of  that  shuttle,  then  we  may  generate  the  new  demand 
rate  by  considering  the  demand  spike  to  occur  at  the  same  relative  posi¬ 
tion  in  the  cycle  (a  proportional  transformation)  and  then  increasing  the 
demand  rate  on  each  day  by  the  demand  rate  of  a  day  in  the  original  cycle 
on  which  there  was  no  spike.  This  will  result  in  the  cycle  illustrated 
in  Figure  1.4. 

rate  Xg(t)  —  —  — 


Figure  1.4  Compressed  demand  in  a  one-shuttle  system. 

This  pattern  has  some  important  similarities  to  the  launch  cycle  shown  in 
Figure  1.3.  Both  cycles  experience  an  equal  number  of  expected  backorders 
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per  cycle  of  length  L/2,  demand  spikes  of  the  same  magnitude  at  equally 
spaced  intervals,  and  interlaunch  cycles  which  are  probabilistic  replicas 
of  one  another.  In  fact,  the  only  important  difference  between  the  two 
demand  patterns  is  the  change  in  the  position  of  the  demand  spike  within 
the  shortened  cycle.  Each  time  the  cycle  is  compressed,  the  maintenance 
day  will  move  closer  to  the  beginning  of  the  cycle  for  the  single  shuttle, 
whereas  this  will  not  be  the  case  for  many  shuttles  with  overlapping 
launch  cycles.  Since  the  weight  of  a  backorder  near  the  end  of  a  cycle 
is  great,  items  with  later  maintenance  dates  will  be  relatively  more 
critical  than  items  maintained  earlier  in  the  cycle.  As  the  cycles 
become  shorter  in  the  one  shuttle  case,  certain  items  will  consistently 
become  more  important  than  others.  This  is  not  the  case  if  there  is 
more  than  one  shuttle  operating:  in  that  scenario  an  item's  importance 
will  depend  upon  how  many  shuttles  there  are. 

The  multi -shuttle  case  is  the  more  realistic  of  the  two,  but  one 
shuttle  with  an  increasing  launch  rate  is  simpler  to  model,  and  was 
sufficient  for  our  analysis  when  the  actual  day  of  the  demand  spike  was 
unavailable.  It  also  represents  a  worst  case  analysis,  with  certain 
items  becoming  more  and  more  important  as  cycles  shortened.  The  results 
of  this  effect  will  be  discussed  in  greater  detail  later.  Lastly,  as  we 
will  demonstrate,  the  effect  of  performing  maintenance  on  one  day  in  the 
cycle  diminishes  as  the  launch  cycle  shortens.  In  other  words,  the  dis¬ 
tribution  begins  to  approximate  a  stationary  Poisson  distribution,  render¬ 
ing  the  position  of  the  demand  spike  extraneous  to  the  spares  determination 
problem. 

The  above  discussion  of  pre-launch  maintenance  operations  may  con¬ 
veniently  be  translated  into  a  single  location  model  concerned  with  the 
inter-launch  maintenance  for  a  single  shuttle.  If  we  allow  the  shuttle's 


launch  rate  to  vary,  we  simulate  the  effect  of  changing  the  number  of 
shuttles  when  each  has  a  constant  launch  rate.  We  are  interested  in 
identifying  a  spares  mix  for  serviceable  items  maintained  during  the 
cycle.  These  stock  levels  should  give  the  best  system  performance  at  a 
given  level  of  investment.  We  require  a  method  of  evaluating  the  expected 
performance  of  the  system  which  recognizes  the  special  demand  distribution 
we  have  described  and  which  allows  for  a  changing  backorders  penalty  over 
time. 

Review  of  the  Literature 

This  problem  is  similar  to  one  encountered  in  setting  spares  levels 
for  a  variety  of  operations  involving  high  cost,  low  demand  items.  Feeney 
and  Sherbrooke  [5]  developed  models  for  systems  experiencing  compound 
Poisson  demand  and  developed  measures  of  supply  performance  that  could  be 
used  to  minimize  backorders  for  a  single  item  and  arbitrary  resupply 
distribution.  These  same  authors  later  optimized  supply  system  per¬ 
formance  under  budget  limitations  through  use  of  Lagrange  multipliers, 
still  for  compound  Poisson  demand  [4].  They  obtained  results  when  system 
performance  was  in  terms  of  proportion  of  demands  filled  by  on-hand  in¬ 
ventory. 

The  model  developed  by  NASA  for  its  own  use  in  spares  provisioning 
had  a  different  objective.  NASA's  model  computed  stock  levels  so  that 
each  item  had  at  least  a  .95  probability  of  being  filled  by  on-hand 
inventory.  It  assumed  stationary  Poisson  demand  rates.  The  NASA  model 
contained  no  analysis  of  an  item's  contribution  to  system  performance 
relative  to  its  cost  [l 2  J - 


Mitchell  [8]  made  changes  in  the  NASA  computer  program  so  that  it 
implemented  marginal  analysis  and  then  computed  the  overall  probability 
that  demands  would  be  met  by  the  supply  system  using  on-hand  inventory. 

By  examining  the  results  for  a  group  of  items  selected  from  the  shuttle's 
avionics  subsystems,  he  identified  those  items  that  accounted  for  the  top 
80%  of  the  group's  spares  costs.  Mitchell's  model  did  not  consider  the 
nonstationarity  we  have  discussed  in  demand  patterns,  nor  did  it  account 
for  an  increased  backorder  penalty  near  the  end  of  the  cycle.  Neither 
the  NASA  model  nor  Mitchell's  model  recognized  that  the  distribution  of 
item  failures  is  often  independent  of  the  number  of  flying  hours  in  the 
previous  mission. 

The  launch  cycle  described  here  was  fully  developed  by  Mucks tadt  [9  ] 
He  defined  criteria  by  which  system  performance  could  be  evaluated  and 
described  algorithms  by  which  these  objective  functions  could  be  optimized 
These  criteria  include,  among  others,  the  total  weighted  probability 
that  demands  are  filled  with  on-hand  inventory,  the  total  weighted 
expected  number  of  backorders,  and  the  total  weighted  expected  number  of 
backorder  days.  These  algorithms  employed  marginal  analysis  or  Lagrange 
multiplier  techniques  to  maximize  system  performance  relative  to  the  given 
objective  function  over  a  range  of  budgets.  The  models  had  not  yet  been 
implemented  or  compared  with  the  methods  of  setting  spares  levels  des¬ 
cribed  above,  and  so  our  objective  here  is  to  determine  which  spares 
stocking  model  one  should  use  in  a  given  situation. 

We  now  have  a  specific  set  of  assumptions  with  which  we  can  model  the 
spares  stocking  problem.  The  models  developed  by  Muckstadt  seem  to  most 
nearly  approximate  shuttle  operations  as  we  have  described  them.  In 
Chapter  2  we  will  continue  to  develop  our  statement  of  the  problem  and 


introduce  the  algorithms  proposed  by  Muckstadt  for  determining  an  optimal 
spares  mix. 


CHAPTER  2 


Given  a  set  S  =  (s-j .Sg,. . . .s^) ,  where  s..  represents  the  number  of 
spares  of  type  i  in  the  spares  mix  S,  we  are  first  faced  with  the  question 
of  how  to  evaluate  its  expected  performance. 

The  model  in  use  by  NASA  gives  the  following  rule:  choose  a  spares 
mix  which  sets  the  minimum  probability  that  any  one  item's  demand  will 
be  met  with  on-hand  inventory  greater  than  or  equal  to  a  constant,  denoted 
by  PCNST.  That  is,  choose  S  as  follows: 

si 

l  P{R.  =  x}  >  PCNST,  all  s.  e  S 

x=0  1 

where  Ri  equals  the  demand  over  resupply  time  for  an  item,  and  is 
assumed  to  be  a  Poisson-distributed  random  variable.  PCNST  may  presumably 
be  varied  to  produce  different  sets  S.  Note  that  we  have  no  way  of 
evaluating  system  performance,  except  by  stating  the  N  values  of 

f  P{R-  =  x>. 
x=0 

If  we  also  consider  that  each  item  i  has  a  cost  c^  associated  with 
it,  we  may  use  a  Lagrange  multiplier  technique  which  begins  with  the 
product  of  each  item's  probability  of  sufficiency  and  attempts  to 
maximize  this  quantity  for  a  given  budget  level.  It  chooses  a  spares  mix 
S  which  may  be  evaluated  by  using  the  following  function: 

n  si 

P0Sc  =  n  [  l  P{R.  =  x}],  all  s.  e  S 
s  i=l  x=0  1  1 

where  P0Ss  is  called  the  system  probability  of  sufficiency. 
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In  this  section  we  will  present  three  different  objective  functions, 
including  total  weighted  probability  of  sufficiency,  total  expected 
weighted  backorders,  and  total  expected  weighted  backorder-days.  Each 
of  the  three  models  will  then  treat  only  one  of  these  measures  of  per¬ 
formance,  and  hence  we  expect  that  they  will  result  in  different  spares 
mixes  S. 

For  example,  we  may  set  stock  levels  for  items  A  and  B  in  Table  2.1 
using  the  NASA  model,  the  Lagrange  multiplier  technique,  and  the  total 
weighted  probability  of  sufficiency  model,  which  we  will  derive  later. 

Without  detailing  the  actual  computations,  we  present  the  resulting  spares 
mixes  S  =  (s^,Sg)  in  Table  2.2  when  the  spares  budget  is  $1,300,000  and 
the  shuttles  are  launched  every  35  days.  We  can  see  that  the  NASA  model 
makes  a  different  choice  of  S  than  the  other  two.  In  this  case  the 
Lagrangian  technique  and  the  weighted  P0S  model  give  the  same  stock  levels, 

but  as  can  be  seen  from  Table  2.2,  they  give  different  objective  function 
values.  For  example,  the  Lagrange  multiplier  technique  evaluates  its  ob¬ 
jective  function,  P0Ss,  at  .621.  In  general,  we  will  use  the  objective  func¬ 
tions  presented  in  this  chapter  to  evaluate  the  spares  mixes  given  by  all 
other  models. 

Table  2.2  illustrates  how  different  assumptions  about  the  spares 
system  can  lead  to  the  use  of  a  variety  of  models  which  can  in  turn  give 
different  values  of  s-j ,S2». .  .  ,sN  to  the  decision  maker.  We  now  present 
those  assumptions  which  were  used  to  derive  the  models  developed  here, 
and  explain  why  they  are  reasonable  in  the  context  of  the  system  des¬ 
cribed  in  Chapter  1.  We  will  also  discuss  the  formulations  of  the  three 
models,  the  algorithms  used  to  solve  them,  and  their  computer  implementa¬ 
tions. 

We  will  present  three  alternative  mathematical  formulations  of  the 
serviceable  spares  mix  problem.  They  relate  to  the  first,  second  and 
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Table  2.1  Item  Descriptions 


Item 

Repair  Time 

Cost 

Daily  Demand  Rate 

Maintenance  Day 

2  B 

60  days 

$415,000 

.0364 

34 

[9HH1 

60  days 

$230,000 

.0075 

3 

Table  2.2  Spares  Mixes 


Item 

NASA 

Lagrange 

Multiplier 

Weighted 

P0S 

A 

3 

2 

2 

B 

0 

2 

2 

Objective 

function 

value 

(.638,  .628) 

.621 

.529 

fourth  of  five  models  discussed  by  Muckstadt,  and  we  will  keep  our 
terminology  consistent  with  his  by  denoting  them  as  models  A,  B  and  D, 
respectively J  Although  the  problems  were  initially  formulated  using 
both  continuous  and  discrete  objective  functions,  we  will  state  them 
as  discrete  models  with  the  basic  time  period  of  one  day.  Model  A  seeks 
to  maximize  the  probability,  weighted  for  each  day  in  a  vehicle's  pre-launch 
cycle,  that  no  item  will  experience  more  demands  during  its  resupply 
time  than  there  are  spare  parts.  This  is  a  weighted  measure  of  the  total 

1  Muckstadt  also  introduced  a  model  to  minimize  the  wieghted  sum  of 
shortage  incidents  (model  C)  and  a  model  to  minimize  the  maximum 
expected  weighted  delay  days  for  any  item  (model  E).  Detailed  dis¬ 
cussions  of  both  appear  in  reference  [9  ]. 


i 


20 


! 

I 


probability  of  sufficiency.  Model  B's  objective  is  to  minimize  the 
expected  weighted  number  of  unit  backorders  for  all  items  and  each  day 
during  a  vehicle's  maintenance  cycle.  Lastly,  in  model  D  we  minimize 
the  expected  weighted  number  of  days  a  unit  is  backordered,  totaled  over 
all  items  and  all  days  in  the  vehicle's  maintenance  cycle. 

In  all  three  cases,  the  launch  cycle  may  be  assumed  to  last  exactly 
L  days,  with  w(j)  corresponding  to  the  weight  applied  to  the  objective 
function  for  the  given  model  on  day  j  (probability  of  sufficiency, 
expected  backorders,  or  expected  backorder  days).  All  models  use  the 
constant  C  to  represent  the  amount  of  investment  available  for  serviceable 
item  spares. 

The  following  information  is  required  for  each  of  the  N  items  that 
compete  for  the  limited  spares  budget: 

ni  =  number  of  identical  units  of  item  i  aboard  one  shuttle, 
c.j  =  total  procurement  cost  for  one  unit  of  item  i, 

T.  =  a  constant  resupply  time  for  item  i,  including  trans¬ 
portation  time, 

v.(t)  =  failure  rate  (failures  per  day)  for  each  unit  of  item  i, 
assuming  one  mission  has  occurred  since  its  last  main¬ 
tenance  period,  evaluated  at  time  t, 
m.j  =  day  on  which  maintenance  crews  prepare  all  units  of  item 
i  for  launch. 

In  all  of  models  A,  B,  and  D,  the  following  assumptions  apply: 

1.  The  group  of  n^  identical  units  of  item  type  i  may  be  considered 
to  be  a  single  item  whose  failures  have  a  nonstationary  Poisson 
distribution  with  parameter  A^(t)  =  n^v.(t). 
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2.  Successive  cycles  of  varying  length  L  are  treated;  there  is  no 
overlap  between  cycles. 

3.  All  items  are  due  on  the  day  they  fail.  This  assumption  serves 
only  to  simplify  computations.  It  may  easily  be  relaxed  within 
the  framework  of  the  algorithms  we  present. 

4.  The  failure  distribution  is  not  significantly  different  for 
items  on  board  the  shuttle  when  it  lands  versus  items  replaced 
during  the  launch  cycle. 

5.  There  is  no  slack  time  between  a  launch  and  the  beginning  of  the 
next  maintenance  cycle  (i.e.,  flying  time  is  assumed  to  be  zero). 

6.  Since  the  days  on  which  maintenance  is  performed  on  the  items  are 
unknown,  the  values  of  m^ ,  which  are  identical  throughout  all 

the  models  for  a  given  cycle  length  L,  were  sampled  from  uniformly 
distributed  random  variables  over  the  range  (0,L). 

Assumption  1  follows  from  the  fact  that  the  ni  units  of  item  type  i 
have  independent  identically  distributed  nonstationary  Poisson  failure 
distributions.  So  for  every  day  j  during  the  cycle,  the  sum  of  the  n^ 
demand  rates  gives  the  demand  rate  for  all  units  of  type  i.  We  showed  in 
Chapter  1  that  there  is  a  similarity  between  a  system  having  decreasing 
cycle  length  for  one  shuttle  and  one  having  an  increasing  number  of 
shuttles  operating  simultaneously,  their  demand  distributions  effectively 
superimposed  to  give  an  overall  demand  distribution.  Since  we  do  not 
actually  know  the  day  of  scheduled  maintenance  for  each  item,  we  may  choose 
either  of  the  two  systems  as  a  model,  and  we  will  take  the  first.  This 
has  the  result  that  given  a  cycle  length  of  L  days  and  maintenance  day  m^ 
for  item  i,  a  second  cycle  of  length  aL  will  schedule  item  i  for  mainten¬ 
ance  at  around  time  am.  (this  is  inexact  due  to  the  fact  that  am^  may  not 


be  an  integer). 

The  backorder  penalty  will  begin  immediately  when  an  item  is  back¬ 
ordered.  This  reflects  the  time-critical  nature  of  shuttle  operations 
as  well  as  our  lack  of  complete  information  about  the  system  and,  conse¬ 
quently,  our  preference  for  conservative  analysis.  Another  area  where 
information  is  lacking  is  that  of  failure  rates  for  newly  installed 
items.  We  have  noted  that  our  information  on  failures  is  tied  to  the 
execution  of  one  mission  since  the  previous  overhaul,  and  so  the  dis¬ 
tribution  of  more  than  one  failure  in  a  given  vehicle  has  a  parameter 
slightly  lower  than  X^(t).  Table  2.3  shows  the  probability  for  various 
numbers  of  failures  over  a  60-day  lead  time  for  an  item  having  a  failure 
rate  over  the  lead  time  of  0.3.  This  corresponds  to  the  highest  shuttle 
activity  rate  and  the  highest  item  failure  rate  n.v.j  contained  in  the 
data. 


Table  2.3  Probabilities  of  f  Failures 


f 

"Used  Part” 

"New  Part" 

0 

.7408 

.8607 

1 

.2222 

.1291 

2 

.0333 

.0097 

3 

.0033 

.0004 
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Next,  suppose  that  an  item  fails  and  is  replaced  with  a  part  having 
a  resupply  time  failure  rate  of  .15.  The  individual  terms  for  the 
probability  that  0,  1,  2,  or  3  of  these  new  items  are  demanded  over  the 
entire  60-day  lead  time  are  also  given  in  Table  2.3.  But  the  "new"  part 
will  only  be  "new"  for  the  rest  of  the  cycle  during  which  it  is  installed; 
since  the  cycle  length  is  8  days,  this  corresponds  to  a  small  fraction  of 
the  total  lead  time.  The  analysis  shows  that  the  probabilities  for  "new" 
parts  are  different  from  those  for  "old"  parts,  but  that  as  the  number  of 
failures  increases  the  magnitude  of  the  difference  drops  off  sharply.  Thus 
we  will  assume  that  the  expected  failure  rate  for  a  "new"  item  is  little 
different  from  that  of  a  "used"  item. 

Assumption  5  states  that  slack  time  will  be  brief  between  the  launch 
and  the  next  maintenance  cycle;  if  an  actual  one-shuttle  facility  were 
under  study  we  would  lengthen  the  launch  cycle  by  the  length  of  an 
average  flight  and  allow  zero  expected  demand  on  those  days.  The  last 
assumption  specifies  our  method  of  randomly  assigning  maintenance  days. 

We  emphasize  that  the  maintenance  schedules  are  identical  between  models 
for  each  cycle  length  and  that  if  L'  =  aL,  then  m^  =  am.,  for  all  i. 

For  each  of  the  items  we  consider,  the  models’  output  consists  of 
a  recommended  integer  value  for  the  spares  level  s.,  where  there  are  N 
item  types  in  the  system.  If  each  spare  for  repairable  item  i  costs  c^ 
dollars,  and  the  investment  limit  is  C,  then  we  have  the  following 
constraints: 

N 

j,  Vi  i c 

si  e  {0,1 ,. . . }  i  =  1 ,. . .  ,N. 


We  are  not  actually  given  the  nonstationary  distribution  for  (t), 
but  are  instead  given  the  failure  (removal)  rate  per  flight,  r. . 
Arbitrarily  choosing  to  place  all  of  the  demand  on  the  item's  maintenance 
day  m.  leads  to  the  following  definition  of  v^t): 


vi(t)  = 


fo 


t  f  m. 


ri  t  =  v 


In  our  discussion  of  model  D  we  will  explain  the  variation  of  this  dis¬ 
tribution  used  in  that  model.  Once  v^t)  is  known,  it  is  easy  to  compute 
X.j (t)  using  assumption  1. 

All  three  models  also  involve  an  expression  for  demands  over  a 
resupply  time  beginning  T.-l  days  before  day  j.  The  lead  time  failure 
rate  A(i,j)  was  introduced  in  Chapter  1.  Recall  that: 


A(i,j) 


l  , 

t=j-T.+l 


Ai(j) 


is  an  approximation  of  this  failure  rate.  If  we  then  wish  to  calculate 
p[Rj (j)-k]>  the  probability  that  there  were  k  demands  for  item  type  i  over 
the  interval  [j-T.-+l,j],  we  may  use  the  following  equation: 


P[Ri(j)=k3 


k! 


Again,  implicit  in  this  equation  is  the  assumption  that  all  items  of  type  i 
("old"  and  "new")  have  indistinguishable  failure  distributions.  It  further 
assumes  the  following: 


Model  A 


We  might  wish  to  maximize  the  sum  of  weighted  probabilities  of  zero 
backorders  over  all  days  in  the  cycle.  The  probability  that  an  item 
experiences  no  backorders  is  also  known  as  probability  of  sufficiency. 
This  objective  function  is: 

L  N 

A(s, ,.. .  ,sj  =  l  w(j)  n  P[R,(j)  <  s.]. 

1  N  j=l  i=l  1  “  1 


System  probability  of  sufficiency  is  a  generally  accepted  measure  of 
supply  system  performance,  and  the  above  expression  extends  its  definition 
to  cases  of  nonstationary  demand  over  unevenly  weighted  cycles. 

A  solution  procedure  for  this  problem  was  discussed  by  Muckstadt  [9  ] 
As  the  objective  function  is  nonseparable  when  there  is  more  than  one 
day  j  with  nonzero  weight  w(j),  we  cannot  employ  a  Lagrange  multiplier 
technique.  A  workable  method  first  computes  an  initial  solution  which 
gives  a  very  low  investment  level  and  a  low  value  of  A(s-| ,. . .  ,sN).  The 
algorithm  sequentially  selects  the  spares  whose  contributions  to  the 
objective  function  are  greatest  relative  to  their  costs.  The  value  of 
A(s^,...,s^)  is  improved  every  time  a  spare  is  added  to  the  mix.  This 


2 

We  justify  this  notion  because  the  items  we  analyze  do  have  very  low 
failure  rates,  and  because  the  higher  an  item's  failure  rate,  the  higher 
its  stock  level.  For  A(i,j)  =  .3,  a  high  rate  for  our  analysis,  and 
s-j  =  4,  n^  =  2,  P[R-i ( j )  =  7]  =  3.2  x  10-8.  Thus  we  may  avoid  the 
mechanical  complexity  of  using  a  truncated  nonstationary  Poisson  dis¬ 
tribution.  To  be  strictly  correct,  the  probability  of  more  than  s.+n. 
demands  should  be  exactly  zero  for  the  one-shuttle  case.  ' 


ii 
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is  basically  a  marginal  analysis  technique  and  yields  good,  if  not 
optimal,  results  for  many  applications. 

The  problem  is  formally  stated  in  Figure  2.1  and  a  diagram  outlining 
the  solution  procedure  appears  in  Figure  2.2.  For  additional  information, 
see  Muckstadt's  paper  [g]. 

L  N  „ 

Maximize  >  w(j)  II  P[R. (j)  <  s.] 
j=l  i=l  1  ”  1 

N 

subject  to  l  c^s.  £  C 

s^  0  and  integer,  i  =  1 , . . .  ,N. 

Figure  2.1  Problem  Statement  for  Model  A 


Model  B 

In  this  model  we  are  interested  in  minimizing  the  total  expected 
weighted  number  of  backorders  for  all  items  throughout  the  launch  cycle. 

As  before  we  will  be  evaluating  the  objective  function  for  each  day  and 
then  summing  across  days  and  items,  rather  than  trying  to  measure  expected 
backorders  continuously  through  time.  For  a  single  item  i,  the  expected 
number  of  backorders  on  day  j  is  given  by  the  following: 

E[B.(j)]  =  l  (x-s.)  P ( R . ( j )  =  x), 
x>s. 

assuming  that  items  incur  a  penalty  from  the  first  moment  they  are 


i 


Yes 


Figure  2.2  Flowchart 
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backordered,  where  R^(j)  has  the  same  meaning  here  as  in  Model  A. 
The  problem  is  formally  stated  in  Figure  2.3. 


Model  B 

N  L 

Minimize  £  £  w(j)  l  (x-s.)*P[R. (j)=x] 

i=l  j=l  x>si  1  1 


N  . 

subject  to  £  c.si  £  C 


si  6  {0,1,...},  i  -  1 , . . . , N 


Figure  2.3  Problem  Statement  for  Model  B 


Muckstadt  [g]  points  out  that  this  is  a  separable  problem  and  rewrites 
the  objective  function  using  a  Lagrange  multiplier  e.  Thus,  for  each  i, 
the  object  is  to  minimize  the  following  for  a  given  nonnegative  6: 


r  L 

l  w,(j)  l  (x-s.)  P[R.(j)=x]  +  ec.s. 
Lj=l  1  x>si  11  11 


si  e  {0,1,...}. 


Since  F„.  is  convex  in  s..  [9  ],  we  may  minimize  its  value  by  taking  first 
differences,  and  identifying  the  smallest  s^  for  which  adding  an  additional 
spare  will  cause  F.j(s..)  to  decrease.  In  other  words,  we  are  finding: 
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nlrriSj'.FjtSj)  -  F^s.+l)  <  0} 

L 

=  min{s. :  £  w(j)[(  £  (x-sf  )P{R,  ( j  )«x>  +  0c.  sj 
1  j=l  x>s.  11  11 

-  (  l  (x-(s .+l))P[R.(j)=x]  +  Bc^s.  +  9c,)]  <  0) 
x>si+l  1  1  111- 

L  si 

=  min{s, :  l  w(j)(l-  £  P[R, (j )=x])  <  0c- }. 

1  j=l  x=0  1  1 

Examining  the  above  equation,  we  see  that  it  involves  the  probability  of 
one  or  more  backorders:  ^ 

si 

1  -  l  P[Ri(j)=x]. 
x=0  1 

This  quantity  provides  some  insight  for  choosing  6.  Denote  the  maximum 
acceptable  probability  of  a  backorder  by  (1-a).  Then  we  must  choose  0  so 
that  0*ci  will  always  give  an  acceptable  upper  bound  to  the  weighted 
probability  of  a  backorder.  If  k  is  the  costliest  spare,  Oc^  is  the 
greatest  upper  bound  we  will  create,  and  we  must  insure  that  the  following 
holds: 


L 

l  w(j)*(l-a)  =  0c 
j=l 


The  result  is  to  choose  0  so  that 


L 

I  w(jMl-a) 
j=l _ 


3  As  before,  we  note  that  the  probability  of  more  than  n,+s.  backorders  is 
close  to  zero.  When  only  one  shuttle  is  involved,  thenactual  probability 
of  more  than  n^+si  backorders  is  exactly  equal  to  zero. 
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If  we  then  allow  e  to  decrease,  the  upper  bounds  will  all  decrease,  and 
we  will  need  higher  stock  levels  to  bring  down  the  weighted  backorder 
probabilities.  The  outcome  of  decreasing  acceptable  weighted  backorder 
probabilities  is  to  increase  spares  costs.  Muckstadt  [9]  proposed 
computing  spares  levels  for  all  items  at  decreasing  levels  of  6.  This 
will  result  in  steadily  increasing  spares  costs  until  we  reach  the  desired 
budget  level  C.  The  successive  values  of  6  and  the  corresponding  invest¬ 
ment  required  are  related  in  Figure  2.4. 


T 


Cost 


I 

! 


0 


* 


0 


Figure  2.4  Investment  in  Spares  vs.  Multiplier  0. 


The  algorithm  used  to  implement  the  above  ideas  is  shown  in  Figure  2.5. 
Note  that  the  original  objective  function  for  Model  B  is  not  directly 
minimized.  However,  if  any  cm  is  "close  to"  C,  then  the  {s^  (i=l,...,N)} 
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which  relate  to  cm  will  give  a  near-optimal  solution.  By  changing  the 
increments  of  6,  we  can  generate  any  number  of  pairs  (c,0).  The  algorithm 
presented  here  produces  B  such  pairs  by  setting  9^  =  60/2m,  m  =  1,...,B. 


Model  D 

A  third  model  considered  by  Muckstadt  [ 9  ]  involved  minimizing  the 
total  weighted  expected  backorder-days  across  the  cycle.  The  expected 
waiting  time  per  backorder,  which  is  just  the  length  of  time  that  a 
shuttle  is  delayed  on  the  average,  may  be  expressed  as  E[D- (j)]. 

Here  we  assume  that  the  length  of  a  backorder  is  independent  of  the 
number  of  backorders  outstanding.  This  is  in  keeping  with  the  earlier 
assumptions  of  fixed  resupply  time  and  low  expected  numbers  of  backorders. 

We  now  require  an  expression  for  the  length  of  time  which  a  backorder 
lasts,  E[D.j (j )].  The  following  relationship  holds  in  the  case  of 

4 

stationary  Poisson  demands: 


E[D.(j)] 


E[B.(j)] 


1 

T. 


t=j-T.+l 


x,.(t) 


First,  we  note  that  were  we  to  use  this  approximation  in  Model  D,  there 
is  a  possibility  of  having  a  resupply  time  demand  rate  equal  to  zero 
should  the  interval  (j-T^+l,j)  contain  no  demand  spikes.  In  order  to 
avoid  this  situation,  we  may  redefine  the  unit  failure  rate  v..(t)  as 
follows : 

'4 

This  equation  follows  from  W  =  L/X  where  W  is  the  waiting  time  for  all 
backorders,  L  is  the  number  of  backorders,  and  X  is  the  arrival  rate. 
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v^t)  = 


ti'ra. 


t  = 


m. 


Thus  the  resupply  time  failure  rate  is  less  variable  than  in  the  case 
when  all  demand  is  concentrated  on  one  day  in  each  cycle.  We  continue 
to  calculate  this  value  as  the  sum  of  the  daily  failure  rates  over  the 
resupply  time. 

A  second  important  observation  regarding  this  approximation  is  that 
it  is  independent  of  the  stock  level  s^  and  so  may  be  included  in  the 
constant  term  w(j)  used  in  Model  B.  We  define  w^j),  the  weighted  number 
of  days  that  a  backorder  on  day  j  of  item  i  will  wait,  as  follows: 


wi  ( j ) 


where  w(j)  has  the  same  meaning  here  as  it  did  in  Model  B.  Using  w^(j) 
instead  of  w(j)  in  Model  B  will  give  stock  levels  which  minimize  the 
expected  weighted  backorder-days  for  all  items. 

An  important  change  to  Model  B  will  be  required  when  appropriate 
values  for  8  are  sought,  however.  Since  the  sum  of  the  weights  may  be 
different  for  each  item,  9m  will  now  be  given  by  the  following: 


9„  =  min 
m  • 


^O-a)  I  w.(j)_, 
_ _ 


ci 


We  noted  earlier  that  the  approximation  discussed  above  is  generally 
used  only  in  cases  where  the  demand  distribution  is  stationary  Poisson. 
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In  order  to  test  the  robustness  of  the  equation  when  \^{t)  varies  over 
time,  a  test  was  devised  using  a  computer  simulation.  The  simulation 
will  be  discussed  in  detail  later,  but  for  now  a  simple  description  of 
the  system  will  suffice. 

Figure  2.6  depicts  the  simulated  system,  consisting  of  two  bases 
with  six  and  eight  shuttles,  and  a  single  depot  with  a  repair  time  of  60 
days.  Lateral  resupply  between  bases  is  not  permitted.  There  is  one 
unit  aboard  each  shuttle  with  a  failure  rate  r..  =  .0166,  and  turnaround 
time  is  50  days.  Four  days  prior  to  each  launch,  the  failure  rate  for 
the  items  in  all  eleven  shuttles  increases  for  one  day.  This  is  a  very 
rapid  activity  rate  which,  as  we  will  later  discuss,  detracts  from  the 
effects  of  nonstationarity.  However,  computer  run  time  is  a  limiting 
factor,  and  so  to  observe  many  backorders  it  is  necessary  to  have  a  high 
activity  rate.  Seven  spares  are  initially  provided  to  each  base,  and 
subsequent  spares  shipments  from  the  depot  are  made  on  the  basis  of 
greatest  need.  5 


X 


1 


.083 


Figure  2.6  System  Used  for  Simulation 


5  A  quantitative  definition  of  base  need  will  be  given  in  Chapter  4  and 
is  a  function  of  the  base's  inventory  position  and  of  the  time  remain¬ 
ing  until  the  launch. 


In  tabulating  the  results  of  the  simulation,  which  ran  for  500 
"days",  we  first  record  the  average  number  of  shuttles  grounded,  L".  Then, 
on  every  day,  we  observe  a  value  of  the  lead  time  demand  rate  multiplied 
by  the  actual  waiting  time  experienced  by  a  backorder  on  that  day  (if 
any).  These  values  yield  some  average  value  over  500  days  of  AW.  The 
comparisons  of  the  values  of  L  and  AW  for  an  increasing  nonstationarity 
factor  p  are  shown  in  Table  2.4  and  are  plotted  in  Figure  2.7.  (p  =  2.0 
indicates  that  the  demand  rate  doubled  once  before  each  launch.) 

Table  2.4  Validation  of  Waiting  Time  Approximation 
U  =  AW  with  Nonstationarity  Factor  p. 


p 

L 

AW 

1.0 

0 

0 

2.0 

0 

0 

2.5 

.022 

.025 

3.0 

.048 

.062 

3.5 

.082 

.106 

4.0 

.126 

.182 

Note  that  as  nonstationarity  increases,  backorders  are  overestimated  to 
a  greater  and  greater  extent,  at  least  in  the  simulation  we  performed. 
We  caution  that  the  simulation  parameters  do  not  approximate  the  system 
we  describe  here  with  any  great  degree  of  accuracy.  The  results  of  the 
simulation  do  suggest,  however,  that  as  nonstationarity  increases,  the 
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approximation  tends  to  overestimate  expected  backorders  in  this  case. 

Sherbrooke  [11]  notes  that  for  simple  Poisson  demand  with  parameter 
X..  and  fixed  resupply  time  T^,  the  probability  distribution  of  waiting 
time  D.j  for  item  i,  is  given  by  the  following  relation: 


P[Di  <  d]  =( 


r  s.-l 

l  e‘x  xm/m! ,  d  <  T. 
m=0  1 


.  d  >  T. 


where  x  =  X.*(T..-d).  In  the  case  of  nonstationary  Poisson  demands,  the 
demand  rate  over  time  period  [j-(T^-d)+l ,j],  would  be  required  for  all 
values  of  d,  between  0  and  T\ ,  for  each  day  j  and  item  i.  Denoting  this 
value  as  A(i,j,d),  we  have  the  following  expression  for  expected  waiting 
time; 


Ti 

EDUJ)]  =  l 
1  d=l 


V1 

o  -  I 

k=0 


il 


-A(i,j,d) 


)(A(1.J.d))1 


k! 


The  expression  for  expected  backorder  days  over  the  cycle, 

N  L 

l.  I  E[D.(j)], 

i=l  J=1 

is  a  separable  relation  and  is  convex  in  (s^ • • ,sm) .  We  now  have  two 
approximations,  and  choose  to  implement  the  approximation  introduced  earlier 
by  Muckstadt. 

Using  the  approximation  for  expected  waiting  time  given  by  Muckstadt, 
we  may  write  an  objective  function  for  Model  D.  The  integer  program  for 
Model  D  appears  in  Figure  2.8. 
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Model  D 


N  L 

Minimize  l  l  vUj)  l  (x-s.  )*P[R.  (j)=x] 
i=l  j=l  1  x>si  1  1 


N 

subject  to  l  c-s.  <  C 
i=l  1  1  ” 


s^  >_  0  and  integer,  i  = 


Figure  2.8  Problem  Statement  for  Model  D 


With  the  exception  that  w(j)  is  now  replaced  with  (j)  and  that 
0m  is  derived  differently  (as  outlined  above),  the  algorithm  for  Model  B 
given  in  Figure  2.5  applies  to  Model  D  as  well. 

The  Computer  Models 

The  three  models  just  described  were  each  implemented  using  FORTRAN 
programs  on  an  IBM  370/168  computer  running  with  VM  operating  system  at 
Cornell  University.  The  programs  each  have  essentially  the  same  structure, 
and  were  modifications  of  programs  written  earlier  by  Cogliano  [2  ].  The 
program  structure  is  shown  in  Figure  2.9.  Element  (7)  in  Figure  2.9 
represents  the  core  of  each  program.  Computer  listings  of  these  three 
subroutines  may  be  found  in  Appendix  A.  NASA  has  provided  failure  and 
cost  data  on  items  in  the  avionics  subsystem  [12].  These  items  are 
reduced  to  those  24  that  NASA  computed  as  being  the  most  expensive  from 
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AND  COMPARISON  WITH 
OTHER  SOURCE’  S  STOCK 
LEVELS'  PERFORMANCE 


Figure  2.9  General  Outline  of  Computer  Models 


a  spares  standpoint.  Together  they  represent  80%  of  the  spares  cost  when 
NASA's  probability  of  sufficiency  model  is  used  and  each  item  is  computed 
to  a  .95  probability  of  sufficiency  (i.e.,  PCNST  =  .95). 

Using  the  program  that  NASA  had  developed  we  generated  stock 
levels  based  on  probability  of  sufficiency  and  compared  them  with  the 
solutions  given  by  Models  A,  B,  and  D.  With  stock  levels  from  the  NASA 
program  as  input,  the  programs  for  models  A,  B,  and  D  computed  their 
expected  performance  using  the  objective  function  for  each  model.  We 
were  thus  able  to  establish  a  common  measure  for  comparison  purposes. 


a 


Circled  numbers  indicate  the  order  of  operation. 
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Similarly,  another  set  of  stock  levels  was  produced  by  a  computer 
program  by  use  of  a  Lagrange  multiplier  technique  to  minimize  backorders 
when  demand  was  assumed  to  be  stationary  and  there  was  no  increased  back¬ 
order  penalty  [1  ].  These  stock  levels  are  also  used  as  input  for 
models  A,  B,  and  D  in  order  to  evaluate  the  performance  of  Lagrangian 
analysis  with  respect  to  the  objective  functions  developed  in  this  section. 

Information  concerning  costs,  failure  rates,  and  resupply  times  are 
not  varied  throughout  the  experimental  procedure.  In  order  to  retain 
comparability  with  the  NASA  model,  we  assume  that  two  shuttles  make 
simultaneous  demands  on  the  supply  system.  (Note  that  this  is  the  worst 
case  of  a  two-base  system  in  which  a  common  depot  is  used  and  the  distance 
between  the  bases  is  ignored,  as  no  two  shuttles  could  be  launched 
simultaneously  from  the  same  launch  complex.)  Shuttle  planners  feel 
that  a  reasonable  range  of  launch  cycle  lengths  would  include  cycles  of 
4,  8,  16,  32,  and  50  days,  and  so  runs  were  performed  for  each  of  these 
activity  levels.  The  last  important  variable  we  have  identified  is  the 
weight  of  the  launch  day  as  compared  with  all  other  days  in  the  cycle. 

To  reflect  a  possible  increased  backorder  penalty  on  this  day,  we  let 
the  weight  be  equal  to  either  one,  corresponding  to  an  equal  weighting, 
or  five,  on  the  last  day  of  every  cycle. 

Another  key  element  for  the  three  models  involves  the  shape  of  the 
failure  rate  function.  If  this  rate  is  nonzero  only  on  the  maintenance 
day  for  item  i  during  a  vehicle's  launch  cycle,  the  resulting  nonstation- 
arity  will  be  as  severe  as  possible.  This  is  why  we  chose  earlier  to 
define  the  failure  rate  function  in  this  manner,  except  for  Model  D, 
which,  for  reasons  discussed  earlier,  must  experience  some  minimal  failure 
rate  on  every  day  during  the  cycle.  We  will  now  discuss  the  impact  of 
this  assumption  in  some  detail. 
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We  have  the  following  values  for  X^t): 


xi(t) 


t  =  m. 


otherwise. 


In  Model  D,  this  equality  is  only  an  approximation.  Therefore,  the  lead 
time  demand  A(i,j)  will  be  some  multiple  of  Lt{.  In  Figure  2.10(a) 
the  lead  time  is  slightly  less  than  L,  and  so  A(i,t-|)  is  equal  to  0. 

At  j  =  t2,  however,  we  have  A(i,t2)  =  L*r. .  Similarly,  for  the  longer 
lead  time  Tj  in  Figure  2.10(b)  A(i,t-|)  =  L*ri  while  A( i , t2 )  =  2(1*^). 
Since  A(i,t-j)  and  A(i,t2)  can  differ  from  one  another  by  at  most  L*r^ 
no  matter  how  t-j,  t2  and  are  chosen,  we  have  the  following  inequality: 

|A(i  ,t, )  -  A(i ,t?)  | 

°i - rtTTTt^J - i’.  k=lor2, 

A(i,tk)  >  0. 


It  is  easy  to  see  that  as  increases,  the  variation  in  lead  time  demand 
represents  a  diminishing  proportion  of  the  constantly  increasing  values 
of  A(i,tk).  As  we  will  show  later,  however,  as  few  as  2  or  3  demand 
spikes  during  the  lead  time  are  enough  to  erase  the  effects  of  nonsta- 
tionarity  in  the  spares  systems  we  examined.  The  argument  may  be  for¬ 
malized  as  follows:  we  recognize  that  T.  =  nL  +  At  where  At  <  L  and  n  >  0 
and  integer.  It  can  easily  be  shown  that  the  average  value  of  A(i,j) 
over  the  cycle  has  the  following  relationship  to  T. : 

1  L 

LiLA(l,j)  r-(At(n+l)+(L-At)n) 

~  TT  =  nL  +  At  • 

Denoting  this  average  ratio  as  Ui ,  we  see  that: 


we  may  also  show  that  the  following  holds: 
lim  VU^)  =  0. 

n-*» 

Since  the  asymptotic  properties  of  V  will  result  from  either  fixing 
T..  and  decreasing  L  or  fixing  L  and  increasing  T.  >  we  note  that  non- 
stationarity  will  become  less  important  ( V ( U^. )  ■+  0)  if  either  the 
resupply  time  increases  or  the  launch  cycle  becomes  shorter.  U.  can  thus 
be  seen  as  a  measure  of  nonstationarity  with  variability  converging 
to  0  as  n  gets  large.  We  will  see  that  as  becomes  less 
variable  for  all  items  i  as  a  result  of  steadily  decreasing  L,  the 
results  of  Models  A,  B  and  D  more  closely  approximate  those  of  models 
where  stationarity  is  assumed.  Finally,  we  note  that  these  results 
support  our  earlier  claim  that  a  single  shuttle  with  steadily  increasing 
launch  rate  will  require  the  same  spares  support  as  an  increasing  number 
of  shuttles.  This  is  because  both  result  in  the  same  values  for  L  which 
will  in  turn  produce  equal  values  of  U.  for  e^ach  item  i. 

A  sample  table  of  the  input  data  is  given  in  Table  2.5.  Computational 
results  are  discussed  in  Chapter  3.  The  programs'  output  includes,  for 
a  range  of  investment  levels,  the  performance  which  the  algorithm  achieves 
in  terms  of  its  objective  function,  as  well  as  the  spares  mix  it  identi¬ 
fies.  The  output  may  also  include  the  value  relative  to  the  objective 
function  achieved  by  spares  mixes  supplied  externally.  The  sources  of 
these  stock  levels  are  the  NASA  program  and  the  marginal  analysis  program 


discussed  earlier. 
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Table  2.5.  Input  Data 


CHAPTER  3 


Each  model's  performance  is  measured  and  evaluated  by  use  of  three 
analytical  approaches.  The  output  of  each  program  is  first  compared 
with  the  NASA  and  Lagrange  multiplier  techniques  discussed  earlier. 

Thus  we  may  determine  the  effects  of  nonstationarity  and  weighted  back¬ 
orders  with  respect  to  these  two  baselines.  If  Models  A,  B,  and  D 
provide  the  same  levels  of  the  objective  function  for  less  cost  than 
the  other  models,  then  we  have  effectively  exploited  our  assumptions 
about  the  system.  Next,  each  model  is  analyzed  to  determine  its  behavior 
in  terms  of  the  objective  function  as  well  as  the  spares  mix  when  back¬ 
order  weights  increase  and  inter-launch  cycles  lengthen.  In  this  section, 
we  will  refer  to  the  cycle  as  the  time  between  consecutive  launches. 
Finally,  we  contrast  the  spares  mixes  at  comparable  levels  of  investment 
for  each  of  the  three  models.  In  this  manner  we  may  draw  inferences  as 
to  the  consequences  to  the  spares  mix  of  selecting  one  objective  function 
over  another. 

The  computer  codes  for  each  of  models  A,  B  and  D's  main  subroutines 
appear  in  Appendix  A.  The  subroutines  not  shown  are  essentially  the 
same  for  each  model,  and  relate  to  input  and  output  functions.  A  comment 
in  Model  A's  listing  reveals  the  procedure  used  to  determine  lead-time 
demand  rates  unique  to  Model  D,  and  so  the  code  for  the  computations  of 
lead-time  demand  rates  is  omitted  from  listings  of  the  other  two  models. 
In  terms  of  their  own  objective  functions,  the  models  performed  better 
than,  or  as  well  as,  both  the  NASA  and  the  simple  Lagrange  multiplier 
models.  For  cycle  lengths  of  4,  16,  and  50  days,  the  results  appear  in 
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graphical  form  in  Figures  3.1-3.18.  A  table  of  values  for  a  35-day  cycle 
is  given  in  Table  3.1  and  sample  output  for  each  program  appear  in  Tables 
3.2,  3.3  and  3.4. 

Table  3.1  Performance  of  Models  for  a  35-Day  Cycle: 


SPARES 

LAGRANGE 

SPARES 

BUDGET 

POS(S) 

BUDGET 

4699440. 

0.29019 

9938780. 

5552240. 

0.32797 

10833730. 

6606240. 

0.43380 

1  1637780. 

6835740. 

0.44005 

12789580. 

8840880. 

0.56348 

13528530. 

9129880. 

0. 60244 

14722580. 

1 1061780. 

0.76431 

15937830. 

1 1838280. 

0.79978 

16808576. 

14952720. 

0.89015 

17882192. 

18801392. 

19915872. 

20911856. 

21668336. 

22825328. 


SPARES 

LAGRANGE 

SPARES 

3UDGET 

E  ( BO) 

BUDGET 

9205780 

4699440. 

2.5893 

10203780 

5552240. 

2.2029 

10370780 

6606240. 

1.8057 

11387780 

6335740. 

1.7145 

12135780 

3840880. 

1.1364 

12538580 

9129380. 

1.0611 

17881088 

11061780. 

0.  6912 

18797600 

1  1838280. 

0.5762 

27862400 

14952720. 

0.3145 

SPARES 

LAGRANGE 

SPARES 

BUDGET 

E ( BO  DAYS) 

BUDGET 

4699440. 

9. 8615 

3279340. 

5552240. 

9.4755 

10355640. 

6606240. 

7.9368 

13938180. 

6335740. 

7.9070 

16766480. 

3340880. 

6. 4390 

20438992. 

9129330. 

6.2991 

22748496. 

11061780. 

3.7093 

24005792. 

11833280. 

2.8763 

26691936. 

14952720. 

0. 8821 

28949904. 

44155440. 

0.0 

32819488. 

MODEL  A 

SPARES 

NASA 

POS(S) 

BUDGET 

POS(S) 

0.70345 

9205780. 

0.63102 

0.76034 

1C203780. 

0. 68576 

0. 80632 

1  0870780. 

0.72495 

0.85227 

11387780. 

0.77347 

0.87660 

12135780. 

0. 81561 

0.90753 

12538580. 

0.83192 

0.93465 

17881088. 

0.91871 

0.95082 

18797600. 

0.95258 

0.96341 

n  r\r*  n  a 

27862400. 

0.99616 

0.97240 

0.98187 

0.98667 

0.98906 

0.99197 


NASA 

SPARES 

MODEL  B 

E(BO) 

BUDGET 

E(  BO) 

1. 1561 

4074440. 

2.9783 

0.9926 

6835740. 

1.7145 

0.8197 

1 1061780. 

0.6912 

0.6721 

13703080. 

0.3783 

0.5606 

15632220. 

0.2525 

0.5159 

19765072. 

0. 0838 

0. 1808 

22026784. 

0. 0463 

0.  1245 

24727168. 

0. 0246 

0. 0133 

28638688. 

0. 0094 

MODEL  D 

SPARES 

NASA 

E ( BO  DAYS) 

BUDGET 

E(BO  DAYS 

2. 9900 

9205780. 

4.9967 

1 . 9000 

10203780. 

3. 3543 

0.  81  87 

10870780. 

3. 1798 

O.4379 

11387780. 

3. 0145 

0. 1764 

12135780. 

2.0523 

0. 0794 

12538580. 

2. 0062 

0. 0553 

1 783 1 088 . 

0. 4438 

0. 0276 

18797600. 

0. 3556 

0.0161 

0. 0050 

27862400. 

0.  0335 

i 


SUMMARY  OF  TOTAL  ASSETS  AND  THEIR  DISTRIBUTION 


SPARES 

ITEM  STOCK  LEVELS 


PROBABILITY  OF 
SUFFICIENCY 


1  4 

2  3 

3  2 

1 

5  1 

6  1 

7  1 

8  1 

9  1 

1 0  2 

11  1 

12  0 

13  0 

14  1 

15  0 

16  0 

17  1 

18  1 

19  1 

2'  0 

21  1 

22  2 

23  0 

24  1 


0.88443 

0.96333 

0. 94474 

0.96360 

0.90250 

0.95409 

0.97313 

0.97313 

0.90274 

0.98567 

0. 99214 

1 . 00000 

1 . 00000 

0.99214 

0.95131 

0.96900 

0.99810 

0.99270 

0.96926 

0.94857 

0.97814 

0. 98930 

0.98325 

0.99043 


TOTAL  SPARES  INVESTMENT  9072640.00 
SYSTEM  POS  0.64623320 


Table  3.2  Sample  Output  for  Model  A 


SUMMARY  OF  TOTAL  STOCK  LEVELS  AND  EXPECTED  BACKORDERS 


SPARES 

ITEM  STOCK  LEVELS 


TOTAL  EXPECTED 
AVERAGE  SHORTAGES 


1 

8 

0.  10 

2 

8 

0.07 

3 

4 

0.  13 

4 

2 

0.24 

5 

3 

0.  14 

6 

3 

0.32 

7 

2 

0.  15 

8 

2 

0.  15 

9 

3 

0.09 

10 

3 

0.  12 

1 1 

2 

0.  15 

12 

1 

0.  10 

13 

1 

0.10 

14 

2 

0.  15 

15 

1 

0.28 

16 

1 

0.  11 

17 

1 

0.45 

18 

2 

0.02 

19 

2 

0.  19 

20 

1 

0.32 

21 

2 

0.11 

22 

3 

0.07 

23 

1 

0.04 

24 

2 

0.03 

TOTAL  EXPECTED  WEIGHTED  BACKORDERS  0.0463 

TOTAL  SPARES  INVESTMENT  22026784.0 


Table  3.3  Sample  Output  for  Model  B 


SUMMARY  OF  TOTAL  STOCK  LEVELS  AND  EXPECTED  BACKORDER  DAYS 


SPARES 

ITEM  STOCK  LEVELS 


TOTAL  EXPECTED 
BACKORDER  DAYS 


1  4 

2  4 

3  2 

4  1 

5  2 

6  2 

7  1 

8  1 

9  1 

10  2 

1  1  2 

12  1 

13  1 

14  2 

15  1 

16  1 

17  1 

18  1 

19  1 

20  1 

21  1 

22  2 

23  1 

24  1 


0.05 
0.04 
0.06 
0.04 
0.02 
0.  04 
0.09 
0.09 
0.  02 
0.02 
0.01 
0.00 
0.00 
0.01 
0.03 
0.02 
0.  04 
0.05 
0.04 
0.  04 
0.03 
0.02 
0.01 
0.05 


TOTAL  EXPECTED  WEIGHTED  BACKORDER  DAYS  0.818674 

TOTAL  SPARES  INVESTMENT  13938180.0 


Table  3.4  Sample  Output  for  Model  D 


Looking  at  Figure  3.1,  we  note  that  in  a  four-day  cycle.  Model  A 
does  perceptibly  better  than  both  the  simple  Lagrangian  technique  and  the 
NASA  probability  of  sufficiency  model.  For  a  given  budget  level,  the 
improvement  over  the  Lagrangian  method  is  about  1%,  measured  in  terms  of 
weighted  probability  of  sufficiency.  When  the  weight  increases  from  one 
(an  equal  weighting)  to  five,  we  see  in  Figure  3.2  that  little  change  in 
the  models'  relative  performance  is  evident.  Moving  to  a  16-day  cycle 
in  Figure  3.3,  we  see  an  increased  differential  in  performance  level, 
to  about  2%.  This  relfects  the  increased  nonstationarity  effects 
inherent  to  a  longer  cycle.  A  similar  increase  is  evident  in  Figures 
3.5-3. 6  when  the  launch  days  are  increased  to  50  days  apart,  leading  to 
a  performance  increase  of  about  2-4%.  The  curves  for  Model  A  are  derived 
by  executing  the  algorithm  presented  in  Chapter  2  and  stopping  every 
time  $1  million  is  added  to  the  total  spares  cost.  Since  the  data 
points  for  the  NASA  program  and  the  Lagrange  multiplier  technique  are 
also  discrete,  due  to  the  integer  nature  of  the  decision  variables,  it 
would  be  misleading  to  fit  a  curve  through  them.  We  are  limited  to 
comparing  points  which  are  more  or  less  adjacent  to  one  another  and 
drawing  inferences  from  their  costs  and  performances  relative  to  the 
objective  function.  Thus  one  conclusion  we  may  draw  from  the  data  is 
that  as  the  launch  cycle  for  a  single  shuttle  becomes  shorter,  the 
effects  of  the  nonstationarity  we  expressed  in  the  formulation  of  Model  A 
are  less  pronounced. 

Secondly,  it  appears  that  the  backorder  weight  on  the  launch  day 
does  not  affect  the  choice  of  a  spares  mix  as  much  as  does  the  non¬ 
stationary  nature  of  lead-time  demand.  It  should  be  the  case,  however. 
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that  if  only  certain  items  were  given  increased  backorder  weight  and 
if  the  increase  lasted  longer  than  one  day,  those  items  would  necessarily 
be  stocked  to  a  greater  depth  in  a  spares  mix  produced  by  Model  A.  Model  A 
would  then  be  more  responsive  to  weights  than  it  is  under  the  environment 
we  assumed  for  our  experiment. 

The  output  for  Model  B,  assuming  a  four-day  cycle  and  equal 
weighting  of  the  launch  day,  is  shown  in  Figure  3.7.  Here  we  see  a 
difference  of  about  2-4%  between  the  performances  of  Model  B  and  the 
NASA  model,  but  no  apparent  difference  between  it  and  the  Lagrange 
multiplier  method.  The  same  pattern  is  observed  in  Figure  3.8  for  a 
launch  day  weight  of  5,  only  the  performance  of  all  the  methods  has 
dropped.  In  fact.  Model  B  actually  chose  the  same  stock  levels  as  it 
did  for  a  launch  weight  of  1,  suggesting  that  it  is  perhaps  not  possible 
to  increase  system  performance  by  modifying  stock  levels  if  a  model 
considers  only  the  increased  weight  of  the  launch  day  relative  to  cost. 

When  the  cycle  length  increases  to  16  days,  in  Figures  3.9-3.10,  we  again 
observe  that  Model  B  gives  generally  lower  backorders  than  does  the 
NASA  technique,  but  achieves  close  to  the  same  performance  as  a 
Lagrangian  technique.  Even  where  the  cycle  between  the  launches  is  50 
days,  in  Figures  3.11-3.12,  we  fail  to  distinguish  an  improvement  over 
the  simple  Lagrangian  technique  when  we  use  Model  B.  Since  the  Lagrange 
technique  has  a  measure  of  backorders  as  its  objective  function,  we  would 
expect  that  of  the  three  models,  the  results  of  Model  B  should  be  most 
closely  approximated  by  the  Lagrange  multiplier  method.  Therefore,  we 
recognize  that  any  improvement  obtainable  using  Model  B  might  not  be 
evident  among  the  small  number  of  spares  we  examined. 
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Fig.  3.7  Spares  Mix  Comparison  for  Four  Day  Unweighted  Cycle  -  Model  B. 
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Finally,  in  Figures  3.13-3.14,  the  expected  waiting  time  over  the 
cycle  is  shown  for  Model  D,  the  NASA  model,  and  the  Lagrangian  model. 

Even  for  the  four-day  evenly  weighted  launch  cycle,  there  are  signifi¬ 
cant  advantages  to  be  gained  by  choosing  one  spares  mix  over  another. 

As  much  as  a  50 %  reduction  in  expected  waiting  time  for  backorders  will 
accrue  from  using  Model  D  to  set  the  stock  levels.  Most  of  the  reduction 
seems  to  stem  from  a  single  item  which  has  such  a  low  lead  time  demand 
that  few  spares  would  ever  be  in  resupply  and  a  long  delay  would  result 
were  it  ever  backordered.  This  item  is  bypassed  for  low  budget  levels  in 
other  models  because  of  its  low  failure  rate.  As  investment  increases. 
Model  D  and  the  other  models  begin  returning  about  the  same  levels  of 
performance  with  respect  to  expected  weighted  backorder  days.  As  was 
noted  in  the  analyses  of  Models  A  and  B,  the  introduction  of  longer 
cycles  has  the  result  of  both  decreasing  the  required  budget  for  all 
levels  of  performance  and  increasing  the  model's  sensitivity  to  nonsta¬ 
tionary  demand  rates  which  increase  instability  in  expected  lead-time 
demand.  The  result  is  that  in  Figures  3.17-3.18  for  a  50-day  cycle 
there  can  be  as  much  as  a  70%  reduction  in  expected  weighted  backorder- 
days  for  low  budget  levels  when  Model  D  stock  levels  are  employed. 

We  next  turn  to  an  analysis  of  how  the  various  models  allocate 
spares  as  the  budget  limits  increase  while  the  cycle  length  is  held 
constant  at  35  days.  The  launch  day  has  five  times  the  weight  of  other 
days  in  the  cycle.  For  each  of  the  models  we  tried  to  find  three  spares 
mixes  priced  at  $10  million,  $15  million,  and  $20  million,  although  some 
of  the  budgets  were  lower  or  higher  than  others.  Budgets  which  were 
higher  by  $1  million  could  include  from  two  to  five  more  spares  than 
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Fig.  3.13  Spares  Mix  Comparison  for  Four  Day  Unweighted  Cycle  -  Model  D. 
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more  austere  budgets,  and  so  the  fluctuations  in  performance  with  respect 
to  an  objective  function  could  appear  to  be  greater  than  they  actually 
are.  The  budgets  and  performance  levels  are  shown  in  Tables  3.5  and  3.6 
respectively,  while  the  item  stock  levels  for  each  of  the  twelve  spares 
mixes  are  given  in  Table  3.7. 

Close  examination  of  Table  3.7  shows  that  of  the  Models  A,  B,  and 
D,  Model  B  gives  priority  to  the  high-demand,  low-cost  items  (1  and  2) 
sooner  than  Model  A  or  Model  D.  Model  A  tends  to  stock  the  moderate 
demand,  high  cost  items  (7  and  8)  sooner  than  either  of  the  other  two 
models;  but,  as  demand  decreases  slightly  and  cost  rises,  and  lead  time 
shortens  (items  12  and  13)  the  opposite  effect  is  apparent.  The  lead  time 
proves  to  be  an  important  factor  for  Model  D  as  well,  for  as  it  lengthens, 
as  in  the  cases  of  items  15  and  16,  Model  D  invests  in  the  high  cost, 
low  demand  items  more  quickly.  Conversely,  for  a  short  lead  time  (27 
days  for  item  21)  Model  D  buys  less  of  an  inexpensive,  moderately  demanded 
item  than  either  of  the  other  two  models. 

The  overall  trends  in  stocking  policy  seem  to  be  very  similar  for 
Models  A  and  B,  whereas  Model  D  seems  to  stock  at  least  one  item  of  each 
type  by  the  time  it  reaches  budget  level  2.  This  could  be  because  it 
overestimates  backorders  for  the  low  demand  items  under  the  assumption 
of  nonstationary  demand,  as  discussed  earlier.  Another  intuitive  result 
is  that  Model  B  closely  follows  the  stock  levels  set  by  the  Lagrangian 
model,  because  both  are  based  on  a  measure  of  backorders  due  to  lead- 
time  demand.  The  items  used  here  evidently  do  not  vary  significantly 
in  terms  of  expected  numbers  of  backorders,  even  when  the  demand  has  a 
nonstationary  Poisson  distribution.  By  far,  the  most  unusual  spares  mixes 
are  produced  by  Model  D,  which  gives  the  same  spares  mixes  as  Model  A 
in  any  one  budget  only  57%  of  the  time,  and  with  Model  B  only  64%  of 


Table  3.5  Budgets  Used  for  Comparison 


Budget 

Actual  Spares  Mix  Cost,  millions 

Model  A 

Model  B 

Model  D 

NASA 

Lagrange 

1 

10.1 

11.1 

10.3 

10.2 

11.1 

2 

15.4 

15.6 

14.0 

12.5 

15.0 

3 

20.4 

19.8 

20.4 

18.8 

- 

Table  3.6  Performance  Relative  to  Models  A,  B,  and  D  of  Different 
Spares  Mixes 


Spares  Mix 

Computed  By 

Expected  Weighted  Probability  of  Sufficiency 

Budget  1  Budget  2  Budget  3 

Model  A 

NASA 

Lagrange 

.719  .923  .985 
.685  .831  .952 
.764  .890 

Spares  Mix 

Computed  By 

Expected  Weighted  Backorders 

Budget  1  Budget  2  Budget  3 

Model  B 
•NASA 

Lagrange 

.691  .253  .084 
.993  .516  .125 
.691  .315 

Spares  Mix 

Expected  Weighted  Backorder  Days 

yUinpuLcu  Djr 

Budget  1  Budget  2  Budget  3 

0.18 


Model  D 

NASA 

Lagrange 


'0 

5 

0 


0.8; 

2.0 

0.8! 


Table  3.7  Stock  Levels  for  Comparison 


SOURCE 

Model 

A 

Model 

B 

Model 

D 

NASA 

Lagrange 

BUDGET 

1 
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1 
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2 
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2 

3 

1 

2  3 

ITEM 

1 

4 

6 

7 

5 

7 
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4 

5 

4 

5 

6 

5 
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4 

4 

6 
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6 

7 
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4 

5 

4 

5 

6 

5 

6 

3 

2 

3 

4 

2 

3 

3 

2 

2 

3 

2 

3 

3 

2 

3 

4 

1 

2 

2 

1 

1 

2 

1 

1 

2 

1 

1 

2 

1 

1 

5 

1 

2 

3 

2 

2 

3 

1 

2 

2 

1 

2 

3 

2 

2 

6 

2 

2 

3 

2 

3 

3 

1 

2 

3 

2 

2 

3 

2 

3 

7 

1 

1 

2 

1 

1 

2 

1 

1 

2 

1 

1 

2 

1 

1 

8 

1 

2 

2 

1 

1 

2 

1 

1 

2 

1 

1 

2 

1 

1 
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1 

2 

3 

1 

2 

2 

0 

1 

1 

1 

1 

2 

1 

1 

10 

2 

3 

3 

2 

2 

3 

2 

2 

3 

1 

2 

2 

2 

2 

11 

1 

1 

2 

1 

2 

2 

1 

2 

2 

1 

1 

2 

1 

2 

12 

0 

0 

0 

0 

0 

1 

0 

1 

1 

0 

0 

1 

0 

0 

13 

0 

0 

0 

0 

0 

1 

0 

1 

1 

0 

0 

1 

0 

0 

14 

1 

1 

2 

1 

2 

2 

1 

2 

2 

1 

1 

2 

1 

2 

15 

0 

1 

1 

0 

1 

1 

1 

1 

2 

1 

1 

1 

0 

1 

16 

0 

1 

1 

0 

1 

1 

1 

1 

1 

0 

1 

1 

0 

1 

17 

1 

1 

1 

1 

1 

1 

1 

1 

2 

1 

1 

1 

1 

1 

18 

1 

1 

2 

1 

1 

1 

1 

1 

2 

1 

1 

1 

1 

1 

19 

1 

2 

3 

1 

2 

2 

1 

1 

2 

1 

1 

2 

1 

2 

20 

1 

1 

1 

1 

1 

1 

1 

1 

2 

1 

1 

1 

1 

1 

21 

1 

2 

2 

1 

2 

2 

1 

1 

2 

1 

1 

2 

1 

2 

22 

2 

2 

3 

2 

2 

3 

2 

2 

3 

1 

2 

2 

2 

2 

23 
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1 
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1 

1 

1 
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0 

0 

1 

0 

1 

24 
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1 

2 
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1 

1 

1 
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the  time.  By  contrast,  models  A  and  B  returned  the  same  stock  levels 
73$  of  the  time. 

The  choice  of  an  objective  function  and  hence  one  of  the  models  given 
above  has  a  major  impact  on  the  spares  mix  that  one  would  purchase  and 
upon  the  performance  that  the  serviceable  spares  system  can  provide 
relative  to  any  one  measure.  The  final  system  performance  cannot  be 
predicted,  however,  until  the  number  of  shuttles,  their  launch  frequency, 
and  their  maintenance  schedules  are  known.  Vet  another  important  system 
characteristic  from  a  repairable  spares  perspective  will  be  the  presence 
of  two  launch  sites  which  will  have  different  activity  rates  associated 
with  them.  It  will  be  important  to  know  how  the  launch  periods  at  the 
two  bases  correspond  to  one  another,  because  when  the  two  time-dependent 
demand  distributions  are  combined,  the  basic  period  for  the  system  will 
depend  upon  their  interactions.  In  some  cases,  the  number  of  days  to  be 
examined  could  be  as  much  as  one  year;  if,  for  example.  Base  A  were  to 
launch  every  4  weeks  while  Base  B  launched  every  13  weeks.  In  the  next 
chapter,  we  examine  the  two  base  problem  as  it  relates  to  the  location 
problem  for  spares,  as  well  as  its  effect  upon  the  choice  of  a  spares  mix. 


CHAPTER  4 


The  choice  of  an  optimal  mix  of  serviceable  spares  is  not  entirely 
separable  from  that  of  where  the  spares  storage  facilities  should  be 
located.  The  most  obvious  interdependence  is  through  the  increased 
number  of  backorders  (and  backorder-days)  we  would  expect  at  bases  where 
no  spares  are  permanently  stocked.  Alternatively,  if  facilities  exist 
to  store  spare  parts  at  both  bases,  a  model  must  recognize  that  only  one 
of  the  two  bases  will  receive  a  given  spare  after  it  completes  repair 
at  the  depot.  However,  especially  in  the  more  likely  case  where  lateral 
(base  to  base)  resupply  is  allowed,  the  increased  transportation  time 
will  only  be  on  the  order  of  a  day  or  two,  or  only  between  two  and  five 
percent  of  the  transportation  times  we  encountered  in  Models  A,  B,  and  D. 
Even  without  a  lateral  resupply  capability,  this  means  that  if  demand 
is  low,  failures  at  each  base  will  be  rare  enough  to  permit  most  repaired 
units  to  be  returned  to  the  base  at  which  they  originally  failed  and 
so  each  base  may  be  treated  as  a  single  location. 

While  the  choice  of  spares  basing  locations  may  not  overly  influence 
the  optimal  number  of  each  item  to  be  procured,  it  could  still  have  sig¬ 
nificant  impact  on  the  real  time  efficiency  of  the  system.  In  the  event 
that  a  backorder  occurs,  the  expected  waiting  time  could  be  significantly 
altered  by  the  existence  of  a  spares  facility  at  each  base  and/or  lateral 
resupply  capabilities.  The  time-critical  nature  of  shuttle  operations 
will  probably  result  in  the  use  of  a  real  time  inventory  monitoring 
system.  Such  a  system  could  be  used  to  track  the  inventory  position  of 
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those  spares  which  could  adversely  affect  the  launch  date,  and  to  direct 
routine  shipments  from  the  depot  to  the  bases  as  well  as  emergency  or 
precautionary  shipments  from  one  base  to  another.  The  number  and 
duration  of  backorders  could  thus  be  held  to  a  minimum,  but  much  depends 
upon  where  the  spares  facilities  are  located  and  upon  the  conditions 
which  will  permit  or  prevent  a  shipment.  One  interesting  kind  of  spares 
shipment  would  take  place  as  the  launch  date  for  one  of  the  bases 
approaches.  If  no  demands  occur  and  the  launch  takes  place,  the  shipment 
will  normally  be  returned  to  the  sending  base.  Such  a  shipment  will  be 
the  result  of  what  we  will  call  a  launch  critical  event. 

Given  the  many  interrelationships  between  the  supply  system,  the 
launch  rates  at  the  different  bases,  and  the  interbase  shipment  discipline, 
it  is  extremely  difficult  to  develop  analytical  models  which  would  provide 
some  optimal  level  of  fills  or  expected  backorders.  Miller  [7]  presents 
a  model  which  can  be  used  as  a  real  time  decision-making  mechanism  for 
repairable  spares  allocation.  This  model,  which  he  terms  Real  Time  Metric 
(RTM),  compares  the  need  of  each  base  with  the  reluctance  of  the  depot 
whenever  a  supply  event  such  as  an  item  failure  or  a  repair  completion 
takes  place.  The  RTM  generates  quantitative  values  for  base 
need  and  depot  reluctance  as  functions  of  the  state  of  the  system.  If 
depot  reluctance  is  smaller  than  at  least  one  base's  need,  a  spare  is 
shipped  to  the  base  with  the  greatest  need. 

This  model  is  easily  extended  to  incorporate  a  nonstationary  demand 
distribution  and  launch  critical  events.  Miller's  model  has  been  pro¬ 
grammed  in  FORTRAN  on  an  IBM  370/168  at  Cornell  University  as  an  experi¬ 
mental  simulation  by  Cogliano  [2]. 
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Using  this  simulation  as  the  basis  for  our  model  of  the  shuttle's 
repairable  spares  supply  system,  we  restate  some  of  Miller's  assumptions 
and  extend  them  as  follows: 

1.  The  system  consists  of  a  depot  and  two  bases.  Any  units  removed 
from  the  shuttle  due  to  failure  or  suspected  failure  are  sent 

to  the  depot  for  repair. 

2.  A  single  item  is  examined. 

3.  Repair  time  is  fixed.  Transportation  time  from  base  to  depot 
is  considered  part  of  the  repair  time,  and  is  fixed  at  t^  days 
from  the  depot  to  a  base  and  at  t^  days  between  bases.  (These 
assumptions  are  easily  modified.) 

4.  The  item  experiences  a  nonstationary  Poisson  demand  distribu¬ 
tion  with  a  demand  spike  mfa  days  before  the  launch  at  each 
base. 

5.  There  is  a  fixed  number  of  spares. 

6.  The  depot  reluctance  is  zero;  thus  repaired  units  are  shipped 
immediately  to  the  base  with  the  greatest  need. 

7.  Base  need  is  a  function  only  of  its  inventory  position  and  the 
length  of  time  before  its  next  scheduled  launch. 

8.  Backorders  are  assumed  not  to  prolong  the  launch  cycle  so  that 
the  launches  at  a  given  base  are  evenly  spaced,  although 

this  assumption  may  be  relaxed. 

9.  The  launch  rate  at  a  base  is  directly  proportional  to  the  number 


of  shuttles  at  the  base  and  does  not  change  over  the  simulation 
period. 


Assumption  1  ascribes  to  the  depot  those  functions  carried  out  by 
the  contractor  in  the  real  system.  NASA  has  no  control  over  the  contrac¬ 
tor's  operations,  and,  as  discussed  in  Chapter  2,  the  shipment  and  testing 
is  so  time-consuming  that  any  removal  must  be  treated  as  a  demand  on  the 
supply  system.  In  Assumption  2,  we  limit  the  scope  to  one  item,  although 
the  simulation  may  easily  be  expanded.  We  may  allow  resupply  times  and 
travel  times  to  be  exponential,  but  the  fixed  quantities  are  acceptable 
as  well.  Next,  although  demand  is  assumed  to  be  nonstationary,  the 
number  of  shuttles  in  the  system  is  so  high  that  the  effects  of  nonsta- 
tionarity  become  unimportant  as  we  noted  in  Chapter  3.  We  use  high  demand 
rates  because  otherwise  it  would  be  very  expensive  to  run  the  simulation 
long  enough  to  collect  significant  observations  of  backorders. 

Each  of  the  fixed  number  of  spares  is  either  in  stock,  en  route  to  or 
from  a  base,  or  in  depot  repair.  As  Miller  points  out,  the  number  in 
repair  is  beyond  our  direct  control,  and  so  the  model  seeks  to  maximize 
the  number  which  are  in  stock  where  the  need  is  greatest,  and  places 
limits  upon  the  number  of  shipments  and  hence  the  en  route  inventory. 

In  defining  base  need,  we  make  use  of  the  following  definition  of  inven¬ 
tory  position  (IP): 

IP  =  #  on-hand  +  §  on  order  -  #  backordered. 

We  also  would  like  to  increase  a  base's  need  substantially  if  inventory 
position  is  negative  or  if  there  is  an  impending  launch  at  that  base. 

When  a  lateral  resupply  is  considered,  the  need  of  the  supplying  base 
is  computed  and  compared  with  the  need  of  the  potential  recipient. 
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Lastly,  we  make  some  assumptions  about  the  launch  cycle.  It  is  again 
necessary  to  assume  that  the  launch  cycles  are  fixed  in  length,  although 
it  is  not  unreasonable  when  considering  many  shuttles  at  a  base  to  allow 
some  latitude  in  the  launch  date.  The  simulation  is  much  more  flexible 
than  the  analytic  models,  and  could  easily  be  modified  were  this  a  crucial 
assumption.  The  assumption  of  a  fixed  number  of  shuttles  at  a  given  base, 
however,  may  be  very  important.  It  could  be  the  case  that  a  shuttle 
will  take  off  from  one  base  and  land  at  another,  due  to  either  landing 
site  weather  conditions  or  an  emergency  landing.  The  subsequent  change 
in  the  launch  pattern,  although  only  temporary,  could  cause  some 
changes  to  supply  system  performance,  especially  if  it  is  a  rigid  system 
which  stocks  at  only  one  location.  The  computer  simulation  could  be 
used  to  analyze  these  transient  effects;  however,  that  is  not  our 
primary  objective  here. 

In  modeling  the  above  system  with  a  simulation,  several  events  must 
first  be  identified.  The  relevant  events  are  listed  in  Table  4.1. 

In  addition  to  the  actions  which  are  produced  automatically  by  these  events, 
there  are  several  actions  which  may  be  produced  by  some  events,  depending 
upon  the  state  of  the  system;  specifically,  the  relevant  variable  is 
the  base  need.  We  now  describe  in  some  detail  how  that  quantity  is 
calculated. 

Miller  defines  base  need  to  be  a  function  of  the  expected  backorders 
over  the  travel  time  from  the  depot  to  the  base.  He  approximates  the 
discrete  conditional  distribution  of  expected  backorders  given  the  on  hand 
inventory  by  a  normal  distribution.  This  technique  is  given  in  reference 
[6].  Then  the  mean  number  of  backorders  is  given  by  the  following: 


Table  4.1  Events  To  Be  Simulated 


EVENT 

ACTION 

1.  Failure 

1 

! 

i 

i 

-  decrement  on  hand  inventory  (-1) 

-  increment  depot  rapair  inventory  (+1) 

-  schedule  a  repair  completion 

-  compute  new  base  need 

-  schedule  another  failure  event 

I 

12.  Repair 
j  completion 

j 

-  decrement  depot  inventory  (-1) 

-  choose  a  destination  for  the  repaired  unit 

-  increment  en  route  inventory  to  that  base  (+1) 
(conditional  event) 

-  compute  new  base  need  at  destination 

3.  Arrival  of  a 

spare  at  base 

-  decrement  en  route  inventory  to  the  base  (-1) 

-  increment  on  hand  inventory  (+1) 

4.  Launch  critical 

(Begin/End) 

-  recompute  base  need  at  the  base 

-  schedule  next  launch  critical  (Begin/End) 

5.  Maintenance 

-  reschedule  the  next  failure  event 
(reflecting  increased  failure  rate) 

-  schedule  next  maintenance  date 

I 

I 


Ug  =  On  hand  +  En  route  -  Atd 

where  A  is  the  demand  rate  over  the  travel  time  td-  Further,  the  variance 
of  the  number  of  backorders  is  given  by  the  following: 


T 


If  uD  is  positive,  it  represents  expected  net  inventory  with  no  backorders 

D 


81 


Miller  then  makes  a  correction  for  the  fact  that  the  true  distribution 
is  not  continuous,  as  follows: 


ar  +  ao  +  Vt 


~  (S°l  -  uB 


2a 


aB} 

aB) 


yB  >  o 
yB  <  0 


where  =  /Og  is  the  standard  deviation.  The  final  step  is  the  calcu 
lation  of  the  mean  of  the  backorders  distribution,  using  the  following 
equation: 


,  -  1/2(ur2) 

E[B]  •  o6  •  £  .  6 


BMB 


i  “(l/2)w2 

27  e 


doj. 


We  are  now  able  to  define  base  need  in  a  variety  of  circumstances.  If 
an  item  fails  at  a  base,  then  we  are  concerned  with  the  desirability  of 
immediately  initiating  a  resupply  to  that  base  from  the  other  base.  How¬ 
ever,  if  the  sending  base's  need  will  increase  past  that  of  the  receiving 
base's  during  the  time  it  would  take  to  ship  the  part  back  and 
forth,  then  the  lateral  shipment  should  not  take  place.  In  effect,  we 
are  minimizing  the  maximum  base  need  over  a  short  time  horizon.  We 
will  thus  carry  out  a  shipment  from  base  1  to  base  2  at  time  t  if  the 
following  holds: 


Nr  (t)  >  max  [Nr  (t)] 
B2  ~  t,t+2tb  B1 


where  ND  (t)  is  the  need  of  base  i  at  time  t.  The  difference  between  this 
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model  and  Miller's  model  is  that  we  must  now  look  forward  in  time  towards 
fluctuations  in  need  which  take  place  at  predetermined  times.  The  other 
events  for  which  similar  comparisons  are  performed  are  summarized  in 
Table  4.2. 


Table  4.2  Conditional  Events 

for  Shipments 

EVENT 

CONDITION 

ACTION  j 

1.  Failure 
(base  i) 

i 

Nr  (t)  >  max  [Nr  (s)] 

Bi  se(t,t+2tb)  Bi 

-  decrement  on-hand 
inventory,  base  i 

-  increment  en  route 
inventory,  base  i 

-  compute  new  base 
needs 

! 

2.  Repair 
completion 

i 

! 

i 

i 

max  [Nr  (s)] 

se(t,t+tb+td)  i 

>  max  [Nr  (s)] 

se(t,t+tb+td)  T 

j 

-  ship  to  base  i  ; 

i 

3.  Launch 
critical 
(base  i) 

Same  as  for  failure, 
base  i 

! 

i 

Same  as  for  failure,  | 
base  i 

_ 1 

The  overall  computation  of  base  need  is  derived  only  partly  from  the 
expected  number  of  backorders  at  the  base.  Other  important  considerations 
include  the  presence  of  backorders  and  the  proximity  to  launch  at  the  base 
In  addition,  at  some  maximum  stock  level,  IPm  .  the  base  has  no  need  at 
all.  Thus  we  may  define  the  need  at  base  i  at  time  t  as  follows: 
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if  IP  <  0 

NBl(t)  =<E[B]+ny1  If  0  <  IP  <  IPmax 

where  n  is  a  constant  representing  the  penalty  for  backorders  as  the 

launch  nears,  and  y.  is  a  variable  taking  value  1  if  a  base  is  launch 

critical  and  0  otherwise.  Thus  base  need  might  be  represented  graphically 

as  in  Fig.  4.1.  We  note  that  expected  backorders  E[B]  as  well  as  yi 

will  depend  not  only  on  inventory  position  but  also  on  the  time  period 

over  which  ND  (t)  is  evaluated. 

°i 

The  choice  of  the  constant  n  must  be  made  carefully,  for  a  system 
which  ships  back  and  forth  with  abandon  can  be  just  as  inefficient  as 
one  which  ships  only  when  a  backorder  occurs.  The  best  value  for  n  will 
depend  upon  the  activity  rates  at  the  bases,  the  expected  backorders,  the 
nonstationarity  of  the  demand  distribution,  and  the  number  of  spares  in 
the  system.  For  our  purposes  it  is  acceptable  to  choose  n  by  using  trial 
and  error,  but  in  any  general  use  of  the  method,  a  more  precise  formulation 
would  be  required.  It  could  be  that  a  good  value  for  n  is  obtainable 
in  the  same  way  that  Miller  calculates  depot  reluctance:  as  an  expo¬ 
nential  function  which  decreases  in  the  number  of  spares  on  hand  at  the 
location  in  question. 

The  above  model  relates  to  the  case  where  there  are  two  stocking 
locations  for  repairable  spares,  and  there  is  an  unlimited  capacity  for 
lateral  resupply  between  bases.  In  our  analysis  we  consider  three 
separate  cases.  Case  I  allows  shipments  to  either  base  upon  repair 


84 


completion,  but  does  not  allow  lateral  resupply.  Case  II  allows  ship¬ 
ment  to  the  second  base  only  if  its  stock  of  spares  is  low  and  it  is  close 
to  the  launch  date  at  that  base,  so  that  limited  lateral  resupplies 
are  permitted.  Lastly,  Case  III  corresponds  to  the  stiuation  where  two 
bases,  each  with  complete  stocking  facilities,  are  allowed  to  perform 
lateral  resupplies  as  often  as  is  indicated  by  the  comparison  of  base 
needs. 

The  Computer  Implementation 

The  model  presented  here  was  coded  in  FORTRAN  and  run  on  an  IBM 
370/168  at  Cornell  University  with  a  VM  operating  system.  Most  of  the 
simulation  structure  is  drawn  from  similar  work  done  by  Cogliano  [2]. 
Significant  changes  include  the  addition  of  launch  critical  and  mainten¬ 
ance  events,  the  extension  of  the  time  frame  for  evaluating  expected 
backorders,  the  introduction  of  nonstationary  demand  rates,  and  the  change 
in  shipment  discipline  to  incorporate  the  three  cases  detailed  above. 

The  simulation  has  an  event-scheduling  format,  and  is  modularly 
designed  so  that  each  subroutine  relates  to  a  specific  event.  The  key 
launch  critical  subroutine  is  given  in  Appendix  B.  The  typical  program 
input  is  shown  in  Figure  4.2,  with  a  flowcr.art  representation  of  the 
simulation  given  in  Figure  4.3.  The  item  which  was  used  for  analysis  is 
like  item  3  from  the  group  of  items  we  examined  in  Chapter  3.  The  item 
experiences  .300  failures  per  day  on  the  average  in  the  simulation  whereas 
the  average  failure  rate  for  item  3  in  models  A,  B,  and  D  was  .318.  The 
system  is  assumed  to  consist  of  seven  shuttles  at  base  1  and  five  shuttles 
at  base  2.  Since  shuttle  turn-around  time  is  fifty  days,  there  is  a 


! 

I 


Fig.  4.1  Base  Need  as  a  Function  of  Inventory  Position. 

seven  day  interlaunch  cycle  at  base  1  compared  with  a  ten  day  cycle  at 
base  2.  There  is  one  item  on  each  shuttle,  each  with  a  failure  rate  of 
1/40  per  day  if  stationary  demand  is  assumed,  and  a  resupply  time  of 
60  days.  There  are  25  spares  in  the  system,  15  of  which  are  initially 
at  base  1,  and  the  remaining  ten  are  at  base  2.  Depot  to  base  travel 
time  is  1.0  days,  while  base  to  base  travel  time  is  2.0  days.  The 
simulations  were  "warmed  up"  for  a  set  period  of  time  before  a  series  of 
daily  observations  were  begun.  The  output,  which  reports  the  average 


Figure  4.2  Input  for  Simulation 

OEACH  INPUT  FORMAT  IS  EITHER  1014  OR  10F4.0 
ONUMBER  OF  BASES 
2 

NUMBER  OF  ITEM  TYPES 
1 

NUMBER  OF  TYPE-1,  TYPE-2,  ...  UNITS  AT  BASE  1 
7 

NUMBER  OF  TYPE-1,  TYPE-2,  ...  SPARES  AT  BASE  1 
15 

HUMBER  OF  TYPE-1,  TYPE-2,  ...  UNITS  AT  BASE  2 

5 

NUMBER  OF  TYPE-1,  TYPE-2,  ...  SPARES  AT  BASE  2 
10 

NUMBER  OF  TYPE-1,  TYPE-2,  ...  SPARES  AT  THE  DEPOT 
0 

AVERAGE  TIME-TO-FAILURE  FOR  TYPE-1,  TYPE-2,  ...  ITEMS 
40.000 

AVERAGE  BASE  REPAIR  TIME  FOR  TYPE-1,  TYPE-2,...  ITEMS 

0.0 

AVERAGE  DEPOT  REPAIR  TIME  FOR  THESE  SAME  ITEMS 
60.000 

PROBABILITY  OF  A  BASE  REPAIR  FOR  TYPE-1,  TYPE -2,  ... 

0.0 

TRAVEL  TIME  FROM  BASE-1,  BASE-2,  ...  TO  DEPOT 

0.0  0.0 

TRAVEL  TIME  FROM  DEPOT  TO  BASE-1,  BASE-2,  ... 

AND  FROM  BASE  1  TO  2,  2  TO  1,  ...,  WITH  LATERAL  RESUPPLY 
1.000  1.000  2.000  2.000 
PROCESSING  TIME  FOR  ORDERS  FROM  BASE-1,  BASE-2,  ...) 

0.0  0.0 

OLENGTH  OF  THE  WARM-UP  PERIOD 
150.000 

TIME  BETWEEN  SUCCESSIVE  OBSERVATIONS 

1.000 

NUMBER  OF  OBSERVATIONS  TO  TAKE  (FORMAT  16) 

500 

ENTER  1  FOR  AN  ESTIMATION  OF  VARIANCES 
0 

ENTER  1  IF  FAILURE  TIMES  ARE  EXPONENTIAL 
1 

ENTER  1  FOR  A  TRACE  OF  EVENTS 
0 

SEEDS  FOR  FAILURE  TIME  AND  PLACE  (2F10.0) 

745623964.  235187469. 

SEEDS  FOR  REPAIR  TIME  AND  PLACE  (2F10.0) 

254768137.  647629632. 

SEED  FOR  INITIAL  CONDITIONS  (F10.0) 

645734621 . 

OBASE  ORDERING  POLICY: 

(1)  ORDER  TO  MATCH  THE  UNIT  THAT  FAILED 

(2)  ORDER  TO  REPLACE  THE  PART  JUST  INSTALLED 

(3)  ORDER  TO  REPLACE  THE  PART  THAT  FAILED 

(0)  DOES  NOT  APPLY;  SHIP  USING  MEED  &  RELUCTANCE 
0 
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(  BEGIN  ) 

,3=- 

INPUT  DATA 


Initialize  number  in 
repair,  completion  times, 
and  time  of  first  failure 


Advance  the  Clock 
to  next  event 


Calculate  Table  of 
Expected  Backorders 


What  N 
event  is 
occurring 


if  statistics  event 


otherwise 


1  BASE 

ITEM 

LAUNCH 

REPAIR 

ITEM 

j  ARRIVAL 

FAILURE 

CRITICAL 

COMPLETION 

MAINTENANCE 

Take  action  based 
on  Tables  4.1  and  4.2 


Collect  statistics 
schedule  new 
statistics  event 


1  Mill  U  \ 

exceeded  / 

?  / 

*r 

Figure  4.3  Flowchart 
for  simulation. 

Print  Output 
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number  of  shuttles  grounded,  is  discussed  in  the  next  section. 

Simulation  Results 

Using  the  hypothetical  item  developed  above,  the  computer  simulation 
is  used  to  test  the  expected  performance  of  supply  systems  under  the  three 
cases  we  identified.  Although  the  system's  configuration  and  shipment  j 

limitations  are  varied,  the  failure  rates  and  random  number  seeds  are 
identical  for  all  runs.  This  results  in  nearly  identical  patterns  of 
failures  during  each  of  the  simulation  trials. 

Our  primary  results  concern  the  case  where  there  is  no  nonstationarity 
in  the  demand  distribution.  As  we  are  dealing  with  an  activity  rate  of 
one  launch  every  4.4  days,  the  results  of  Chapter  3  would  suggest  that 
incorporating  nonstationarity  is  unnecessary.  The  simulation  output  is 
shown  in  Table  4.3  for  the  three  cases  defined  earlier.  We  see  that  the 
best  performance  is  returned  by  a  system  allowing  unlimited  lateral 
resupply.  In  addition,  the  output  would  suggest  that  in  this  case  it  is 
preferable  to  have  both  basing  locations  stocked  with  spares  than  to  stock 
only  one  location.  This  fact  is  demonstrated  by  the  70%  reduction  in 
backorders  which  Case  I  achieves  relative  to  Case  II.  Table  4.4  gives 
information  on  the  performance  of  the  three  decision  rules  with  respect 
to  the  number  of  shuttles  grounded  during  launch  critical  periods.  It 
is  possible  to  achieve  close  to  zero  backorders  if  Case  III  is  implemented. 

Taken  together.  Tables  4.3  and  4.4  suggest  than  an  optimal  policy 
in  all  respects  might  be  the  use  of  Case  III.  We  note,  however,  that 
even  when  the  launch  critical  constant  is  zero,  we  will  make  close  to 
100  shipments  of  our  hypothetical  item  in  a  500-day  time  period.  Thus 
Case  I  might  be  a  better  choice  if  the  cost  of  a  lateral  resupply  capa¬ 
bility  is  high. 


Table  4.3  Average  Number  of  Shuttles  Delayed 


Launch  critical 
constant,  n 

I 

Cases 

II 

III 

0 

.066 

.222 

.008 

.01 

.066 

.222 

.014 

.05 

.066 

.222 

.012 

Table  4.4  Average  Number  of  Shuttles  Delayed  During  Launch  Critical  Period 


Launch  critical 
constant,  n 

I 

Cases 

II 

III 

0 

.076 

.01 

wfm 

.208 

.05 

■Hfli 

.208 

■ 

Table  4.5  Average  Number  of  Shuttles  Delayed  with  Nonstationary  Demand 


Launch  critical 

Cases 

constant,  n 

I 

II 

III 

0 

.400 

.500 

.398 

.05 

.400 

.386 

Table  4.6  Average  Number  of  Shuttles  Delayed  During  Launch  Critical 
Period  with  Nonstationary  Demand 


90 


Nonstationarity  in  the  lead  time  demand  rates  may  easily  be  inte¬ 
grated  into  the  framework  of  Miller’s  RTM.  Since  the  simulation  kept 
a  current  value  of  the  next  failure  time,  the  demands  in  a  given  day 
could  effectively  be  doubled  by  simply  combining  the  expected  number  of 
demands  for  two  days  to  produce  the  demand  for  a  spike  day.  We  chose 
to  do  this  on  the  day  the  shuttles  become  launch  critical,  but  it  could 
occur  on  an  arbitrary  day  or  days  in  the  launch  cycle,  with  a  failure 
rate  of  1/60  per  day  on  non-peak  days  and  1/30  per  day  on  spike  days,  we  see 
in  Table  4.5  that  Case  III  is  able  to  completely  avoid  backorders  for  n=-05. 
When  we  examine  the  performance  of  the  three  cases  during  the  launch 
critical  period,  in  Table  4.6,  we  note  that  Case  II' s  performance  worsens 
during  the  launch  critical  phase.  This  is  due  to  the  fact  that  we  may 
ship  a  spare  to  base  2  only  rarely  for  lack  of  prepositioned  spares 
facilities  at  that  location.  Thus  many  backorders  at  location  2  will 
last  at  least  a  day  or  two.  In  this  case  we  note  that  increasing  the 
launch  critical  constant  in  Table  4.6  may  remedy  the  situation  for  Case  II. 
It  has  no  effect  on  either  Case  I  or  Case  III.  It  is  interesting  to  note 
that  the  precautionary  shipments  which  we  allow  in  Case  III  are  able  to 
eliminate  all  backorders  during  the  launch  critical  phase.  However,  we 
must  caution  that  these  data  may  only  be  considered  as  preliminary 
results  since  they  are  based  on  single  simulation  runs  of  500  days.  Much 
more  investigation  is  necessary  before  we  can  actually  identify  a  good 
value  for  n  in  all  cases  and  before  we  can  state  that  Case  III  is  always 
superior  to  the  other  cases  with  the  proper  choice  of  n- 

Once  the  desireability  of  lateral  resupply  capability  and  multiple 

location  spares  prepositioning  is  determined  for  all  items,  decisions 

as  to  the  shipment  discipline  for  the  overall  system  may  be  made. 
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This  determination  will  have  to  depend  upon  the  probable  utilization 
of  the  lateral  resupply  system  relative  to  its  cost. 

There  is  thus  no  single  relationship  that  determines  which  basing 
concept  is  best  for  all  spares  in  the  sense  that  it  simultaneously  mini¬ 
mizes  grounded  shuttles  on  the  average  as  well  as  during  launch  critical 
periods.  The  launch  critical  constant  as  well  as  the  shipment  discipline 
must  be  carefully  chosen  using  the  best  information  available  about  each 
item.  Methods  such  as  those  outlined  in  this  chapter  can  be  invaluable 
in  performing  the  necessary  analysis  and  in  implementing  the  results. 
However,  one  might  well  attempt  to  formulate  some  combination  of  Models 
A,  B,  and  D  with  the  simulation  introduced  above.  It  may  be  that  for 
the  price  of  lateral  resupply  capability  we  could  substantially  enrich 
the  spares  mix  and  surpass  the  performance  of  Case  III.  The  objective 
of  such  a  formulation  would  be  to  minimize  the  combination  of  spares 
investment  with  the  expected  ongoing  costs  of  transportation  and  spares 
facilities.  In  addition,  the  flexibility  of  a  multi -base  capability  in 
case  of  a  real  emergency  merits  additional  consideration. 


CONCLUSION 


The  results  of  Models  A,  B,  and  D  point  out  the  potential  differences 
in  spares  mixes  and  performance  levels  which  result  from  the  choice  of 
any  one  model.  They  indicate  that,  at  least  for  the  items  we  considered, 
nonstationary  Poisson  demand  rates  may  not  yield  significantly  different 
spares  mixes  than  stationary  rates.  As  lead  time  demand  variability 
increases,  however,  nonstationarity  could  be  an  important  factor  in 
setting  spares  levels. 

The  simulation  results  indicate  that  a  good  real  time  shipment 
policy  should  respond  to  nonstationary  demand  rates  and  to  launch  critical 
events.  Further,  systems  with  both  lateral  resupply  and  prepositioning 
of  assets  at  both  bases  seem  to  return  the  best  performance.  The  final 
judgments  may  not  be  made,  however,  until  the  cost  of  each  of  these 
capabilities  are  included  in  the  analysis. 
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APPENDIX  A  -  MODELS  A,B,  AND  D 
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STAPES  MX  DETERMINATION  AiGOSITHK  (ft AIN) 
E Y  KATHLEEN  CONLEY,  DECDMEEE  19B1 


c 

c 

c 

c 

c 

c 

c 

c 
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c 

c 


DEECFiif  LICK; 


THIS  IE  THE  /.AIK  SUBROUTINE  FCF  HO  LEI  A.  SIMI¬ 
LAR  DULROUT  IMS  Ea  r  CUT E  MODELS  E  AND  D.  MANY 
EAR  &KEI  E»S  A  ZD  FUNCTIONAL  VALUES  ARE  PASSED 
laF-CU  3h  THIS  SOESGUIISI.  FIRST  THE  SYSTEM  Cl'AS- 
AOl  Zi  LSTICS  AM  Dx.AL  IN  (xiiE'J'l)  ,  ICxiCifxl  -.  r 
Ma'ES  I  I  CM  AN  ZMIiNAl  SOURCE.  IDEE.  LI  i  PEE 
DIMS  AIL  L VALli A" DL  D  Y  MODEL  A  (COM)  ADD  SPALLS 
DIALS  FECK  MCLLL  A  ARE  GENERATED  (AA).  FINALLY, 
IKE  S.GCK  LEVE1S  AND  PELF  CM  ANCE  LEVELS  rOP  ALL 
STAFFS  SIXES  ALL  LiSFLAYED  (OU‘I)  •  VAiLAbl- 
D  E  r  L  K 1 1 1 0  N  S  A  At  AS  F CLICKS: 


boiil  (1) — 1  HE  AVERAGE  DAILY  LEKAliD  fCr  ILLS  1  (INCLUDING 

SEINES- ) 

SJX  (I)  -  -  A  IIZASURE  Cf  FaPECxED  sliGELAGIx  RGB  STATIONARY 
L  i x,  A  K  L  . 

jJCLSI  (1)  --FIXED  RESUPPLY  LIME  FCE  I ‘I  EE  I  (IN'  hOULS) 

LAF.LDA  (I) — M INIS U ft  ACCEPTABLE  STOCK  LEVEL  PCS  1 1 EM  I. 

*.(.!)— TnL  ir,  EIGHT  CP  A  EACKCFDEE  OK  DAI  0. 

EGOS: (I) --ITLX  I  COST  (IN  DOLLARS) 

SiiCKY— A  SErtSURE  01  YC1AL  EXPECTED  SiiCFIAGES. 

EAC  (-.)  — PRCfcAEiLllY  OF  A  E ACKCRDER  Or  llLft  x 

COST  — C  UNBENT  CCS!  Cf  A  SPARES  MIX. 

EOS  U) — FI’ CD  ABILITY  CF  SUIrxCIENCY  -•  0 t-  I  ILL 
ON  THE  LASI  LAY  Cf  THE  CYCLE. 

LUELAY — li.DE  EEIXEE*  ITEM  FAILURE  AND  LATE  DEEDED 
(ASSUME  =  C ) 

LODGE  I  —  MAXIMUM  EDDGE1  LEVEL  FOR  SPAMS  MX  CURRENTLY 
t  L i  i,  j  CCLFUTLD  . 

LSiiCI.L  —  A  XEASOf  E  CF  101 AL  EXPECTED  AVER  AGE  SHORTAGES. 

CExLINo  ON  EUDGa-  1  w .  A (* Y  j ? A , ^ r S  ft  x  a  • 

OPTION  —  SI G K  M  x  0 P  *.I_Trv EE  GCXi'  A—  xsC„  xS  xxxNu  KuLr 
Or,  viHELHxE  MODEL  A  IS  COMPUTING  SPA:  IS  NIXES. 

ICTEDL — TOTAL  FXPiClED  DEMAND  PER  CYCLE  (ASSUMED 
ST  ATICKAP.il  Y) 

::i  I  ■  “  «  x  y  » ,t—  For.  I.  x  x  iv  Ei*  L  E  s  C  i  1  iji,  X  E  ft  1,  „  X  fi  D  F  ii  r.  i,  C  A  I  C  U 


L A  TDD. 

i,  i  x  u  i.  i  —  i.  E  lui.T  OF  F  ACF.CRxLF.  S  ON  L  A  U  K  x  i.  DAYx. 

no  — SU.I  Or  a  A  C"  ?,  CEDE  R  E  h  0  E  A  E  i  L  i  T  L  i_  x-  (NOT  n  P  L  F  r  C  i  .  xi  A  N  0  x 


ft  F  A  S  U  P  F . ) 

SiSiCL  —  POS  (S)  rOP  A  GIVEN  SIMMS  1. 1 X 

NLA  Y  — 1  >,iCi  Till  NUDES?  OF  LAYS  ElT.ELN  LAUKCrixS. 

;  r  —  —  Y  ■.  x  s  x  1  i.  r  it  x  .  i  L  x  x  O  x  E  x.  ft  A f.  x  S x'  L r,  ^  S  x  K  \  x  *  x  D  x  ft  x  Y  C  L  x 

L  x  a  (  .  .  —  CD.  ;  i  x  -  x,  0  1  x  x  V  z  Y  0  *  1 1  i  —  • 

,i .  „ I  .  .  u )  “*  L 1.  1 1.  _i.  x  .  x  x  x  *.  ,(i  ,t  x  r  ri  — .  E  x  0 ...  —  _  _  *  —  0 .  x  A  Y  D  ♦ 

..  —  T  . .. . ; -  -  . .. :  T  v.  Y  /.  I  „  JxLiv  x  x  ii.jj  x  I —  _  Y  x x  .  x  x 


,  NUT--'.  U  ft  MB  01  SiiUI 
(MSUMD  =2) 


.LS  IN  Ex  0  A 


■1  j- 


.  ,u : 


-ii 


v, 


V  “ 


nnnnnnnn  nnn  nnnnnnn 
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SPIKE  (J) — LAKES  OK  VALUE  *1  IF  DEflADD  SPIKE  CL  DAY  J  — 
C 1  H  E  f.  ftlSE  —  C  • 

LS  F —  E  A I  £  Cl  EIDS1  LAUKCii  (ASSUEEL  =  :!EAY/2) 


*************** 


CCEMGli  3DLE  (ICC), SEX  (ICC)  ,E.DCL3L  (IOC)  , LAE  EDA  (100)  , 
*  ft  (100)  ,BCCS1  (100)  , SHOE!, E AC  (100)  , 

*CCS1 ,FCS  (ICC) 

CO!,  .iv  n  u  0  L  u  a  Y  ,  xU  L  u  ET ,  x  S  E  C  3  i  ,1003x0, 

♦GP'IIC*  ,ZC1  DE~  ,  Ell ,W5IGK1 ,EEC  ,3 YSICS  ,  1*EA1  ,  i;s 
CCXSSE  LEX  (ICC) 

0  x’. u  I*  A  L  A  r.  ( _.  C  ,  1v  v  ) 

CCLKOK  KllEXS, NSEU1, SPIKE  ( ICC)  ,I5P 

*****  BESIK  EXECUTION 

CALL  IK PUT 
CALL  CCLt 
C  A 1 1  A  A 
CALL  CCEPC1 

hide:. 

ENL 


SlELCUllNE  CADE E  (II) 
1NIEGEE  II 


IU1S  S  u  B  I-  C  C 1 1  w  L  CAICU1A1E3  1HE  LEAD  USE  ExEAUD  EOF  II  EK 
il  CD  ALL  LAYS  EUE1KG  llix  CYCLE  OE  LENGlh  t«  D  A  Y . 


************************* 


*********** 


c 

CCK:,Ci.  BDCI.  (  ICC)  ,  SEX  ( 100)  ,  DECK  ST  ( 1  J  J  )  ,  1 A  K  3 1)  A  (100)  , 
*  fc  (ICC)  ,ECCSi  (100)  ,SUC£T,£AC(1C3)  , 

*CCSI , ECS  (100) 

cufl.'iuK  DOLL  A  Y  ,  £  U  D  o  ET  ,  1  S  Ei  c  u  1  ,  1  o  C  c  x  o  , 

♦caici  , LOTDKE,  EEL  ,  X  EIGK1  ,  EEC  ,SY  SPSS  ,  LEA  Y  ,KS 
cc*c;.  i  e  x  (i:C) 

CCCHOK  ALAS  (30,  100) 

C C i-j A 0 a  LULLS,  KSiiCl,  SPIKE  (100)  ,131 
ki  ii  1  L  x  S 
C 


ti=? LCE.31  (ii)  /:-u. 
s  1. 1  j.  l  =  i : . 

E 1 0  =  L  u  E  L  A  Y 


(II) 


0=0. 


i  s  k  =  i  s 


c 

C  ****<•  EDO  ID  C  A  EC  ii  I  All  C  l.S  ICi.  LEAD  T 
cc 


nnnnnnnn 
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DO  50  JJ=1,KCAY 
J=J+ 1. 

IF  (ISK.UE. J ) GO  TO  6 
I5L=lSi;*KCAY/2. 
b  IA  =  o-H+ (lEK-u)  *1. 

I  E= J 

7  IF  (1A.GE.C)  GO  TC  8 

IA=IA+KLAi 
IE=I2+NEAY 
GO  TC  7 

o  C  C  0  T X  N  U  E 

A=0. 

f=i  a 

EC  10  f  .  =  J.  A  ,  I F 
i!  t\-i\  I!  ♦  1 

Xr  (FK.G1. KDAY)  KK=1 
A=A+£F1KE  (KK) 

1C  CONTINUE 

AlAfl  (II,  J)  =  A*  (DC?*NCAY)  /2. 

C 

£****************:**************#*•********  *******  ##;***#****«** 

1 

IN  M  0  C  1 1  2  , 

AIAl-i  (II,  J)  =  A*  (DDB*NDAY)/3  +  ICE/5*II 


*****  CO  tit  UIE  KIN  IE  A I  ACCIET8IE  STOCK  IEVEI  (MODEL  A  ONLY) 


110= A I  AM  (11,0) 

IF  (1IJ.G1.SE1WL) GO  TO  45 
5HIKL=1I J 
45  CONTINUE 

50  CONTINUE 

1 AHLBA  (II)=IKT  (SMINL) 

NET L In 
EM 
C 

SuEECUTlNE  INPUT 

C 


c 

C  AN  AFEFEV1ATED  VEFSION  OF  TEE  INPUT  SUEBCUIINE  IS  SHOAN 
C  ELIO  i\ 

C 


C 


iXSDGi.  D JDii  (10  0)  ,SUX  (10C)  ,  FDCF.ST  (100) 
*  >.  (10D)  ,  ECCST  (100)  ,S(iGFT,EAC  ( ICC) 

*EC£T ,:CS  (100) 


*  o  r  ’i  i.'.o. 
C  L  2-  .0  .» 
COM SON 

Cv.  r. ..  j .« 


D  O  w  i"  A  Y 

r  /’  f-  m  *  < 

i  O  i 


f  l  b  L  U  1  i  /  -I  k0  »"i  O  t\  i  |  i  C  O  l'  1  U  f 

/J  I.  x  f  A  uloli  a.  f  X-  £  O  j  b  1  J  1'  W  h 


*  »• 


IE/(KC) 

ill*  ti  ki  (  J  O  f  1  w  w  ) 


:r.u: ,  zi  ilz  ( ICC) 


#  * 


I  Aj.LJA  (11^)  , 


INTEGER  ,  L  ASTI  ,FIli  ST  J,  LAST  J, £  J  AN  (3')) 

C 

C  *****  guAN(l)  IE  int  QUANTITY  CF  ITEM  I  E  E  F  VEHICLE 

C 

C  *****  DEAD  IK  INPUT  DATA 
C 

BEAD  ( 5 , 6  C 1 )  Slltas, KStiUT,CPTION 
ESAU  (5,602)  EUCGE1, NDAY 
HEAL  (5  ,bC€)  L.S, I£E 
READ  (5, 603)  CUELAY, ALFUA 
BEAD  (5,604)  KEEPS , HEIGHT 
601  EC  EM  AT  (2110,  HO.  C) 

601  FORMAT  (FIG.  C,11C) 

60k  FCFV.AC  (10110) 

603  FCivMAT  (2£1C.  2) 

604  FCFMAT  (11U,f1C.2) 

T0CE1G=  BUDGET  +*5CCC0C0. 

EC  610  1=  1 , K ITEMS 

EE  AO  (5 , 6C5)  3CCST  (I)  ,  BECRST  (1 )  ,  ECDB  (I)  ,  QU  AN  (I ) 

605  I  CIX  AT  (31 10. 3, 1 1G) 

C 

c  *****  CCFEICT  ISSUE  ET  FC h  FAILURES  PEP  HOUR  GnOUND  TIME 
c  *****  FLIGHT  TIME  EZF  FLIGHT  IS  96.  G  HOURS 
C 

ESDI,  (I)  =  E£DE  (I)  *96.C*KSHUI*tUAN  (I)  /  (KoAY/2) 

610  CONTINUE 
IOTDtfD=C. Z 
*  T=  V.  £  1G  FIT 
DC  fc  11  J=1  ,  K  E  A  x 
SPIKE  (J)  =C 
V  [J) =1. 

611  CONTINUE 
1£K=1£P 
C 

c  *****  I  NX  EE  ASS  BACKCfDER  WEIGHT  CN  LAUNCH  DaTZS 

C 

DC  G  1 2  L=1,;:s 

IS  K=ISF.+  (1-1)  *KDAY/KS 
W  ( I S  K )  —  r.  1  *  2  .  /  N  S 
61z  CCKIIKUS 

DC  6  17  1=1, KIT  IMS 

totd:x=tctdkd+bddf  (i) 

6  17  CC  m!  I  is  U  E 
C 

r.ETUFt: 

EMC 


s  u  t  r  u  u  1 1  w  I  gc.,p 


******** 


121S  SLJproUI.Ki.  :T<‘D3  IK’  STCCti  j-  ~  V  "•  '•  A  .<  r.  X  .  r  .  •->  J 

OuUECE  itUD  CA1IS  THE  SUBROUTINE  PCI  .CL  El  A  IK  C»Jii 
EVALUATE  I::EM. 
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C 

COKMOK  BDCE  (  100)  ,SEX  (1G0)  ,F,DCFST  (100)  , LAMBDA  (100)  , 
*  fc  (ICC) ,ECCSI  |100) ,SHCfT,£AC  (100) , 

♦COST , FC3  ( 1 00) 

CCiiKUK  DU^DAY  , EUDGET ,ISHCBT ,TCOFIG , 

*C  PT1C  I, ,  TCTEK  0  ,  r  E  I  ,  W  EiG FT ,  EEO  YSPCS  ,  KD  A  Y  ,  KS 
COKfiOK  LEX  (ICC) 

COfc 2CK  ALAS (3G,  ICC) 

CCK301.  KI1LKS,  ASHUT,  SPIKE  (IOC)  ,ISP 

c 


nil’ll  if  ,3) 

3  lCF.\AT(1ii  ,  '  CO EPEE  ISO*  'mill  EASE  CASE') 

W  E 1 7  £  ( t ,  4 ) 

4  rOI -iA'1  ( 1H  /,3X,  'PCS’  ,5X,  ' EEO • , 5X, • 13 VESTKIKT • ) 
3 SAD  ( i.  ,  1 )  ftl  EC 

1  F  C I  ?.  A  7  (13) 

OP71 JN=2 
Ef 2=0. 

DC  10  1  =  1  ,  Nt.EC 

LEAD  (2/2)  CCS  a  ,  (LEX  (J)  ,0=1,  NIT  ESS) 
a  tOisSAT  (115.  C, 3013) 

CALL  A A 
CALL  OUTPUT 
EF1=1. 

»EIlt(t,5)  £Y  SPCS ,1  EC  ,CCS7 
WFILE  (16,20)  COST , S YSPOS 
20  ICEMAT  (IS.  C  ,ffc.5) 

5  FCSflAl (F 10. t, I  10. 6,F  1£. C) 

10  CCMIKUE 

CF7IGK-  1. 

E  E 1 U  li  t< 

IKE 


C 

c 

c 


SUEECUTIKE  OUTPUT 


:  i-  ifc  *  =  =  * 


c 

C  THE  OUTPUT  SUERGUT282  IS  KOI  IMCLUyLD.  1CL  EACH  EUE GET 

C  GEMEATID  CL  EVALUATED,  THIS  SUBKCUTI KE  OUTPUTS  ThE  HE  I. 

C  STOCK  LEVELS,  TEE  VALUE  CF  ECS  (I)  Gl.  THE  LAST  CAT  IN  THE 

C  CYCLE  ICR  iTEn  1,  AMD  TEE  EXFECTED  PEE F CI.MADCE  PCS  (S) 

C  Of  T  :i  E  ASSOCIATED  SPATES  K1X. 

C 


:**■,***************** 


c 


CCr.SC:.  PL  Li.  (ICC)  , S LX  (ICC)  ,r.LCEST  (ICO)  ,  1 A D EDA  (ICO)  , 

*  (100)  >  ECCS7  (IOC)  .SHORT, EAC  (ICC)  , 

*  c  c  s  i , :  c  s  ( i  c  c ) 

^.v. t  . ■  ■  i ■  -  o  L  A  i  ,  LuLv;Ei,iSh<.. IT  ,  , 

' 1  *  1  *  *  t  *  v>  ^  ^  f  L  F  *  f  a  w  J  o « i  x 

CCECCi;  LEX  (ICC) 

CC/EU,.  MAE  (3  C,  ICC) 


100 


CCKdCfc  KIlEMS,NSHt!l,SFIK£(10C)  ,1SF 
CG K MO”  /STATE/  KEAY  (30) 

W  F  IT  E  (6,5)  SYSECS 
FCEfi  AT  ('PCS  (S)  ^  •  ,15. U) 

EEC=0. 


n  n  r>  n  n 


ini 


SUBI.0U3IUF.  10  J  KB IE ME  K1  A1GCF.ITHM  A»  D  EVALUATE 
SFAfES  MIXES  fCF  MOD El  A 

E  Y  K  Alii  LEIN  CCN1EY,  EECEMEEF  19B1 


S  UERO  UTIN  E  A  A 


C 


c 

C  VARIABLES  UNIQUE  1C  THIS  SUEBCD1INE  ABE  DEFINED  BELOV.: 
C 


C  DSHD — SEED  FOB  BAD DOM  NUMLEB  GZNEF.ATCF 

C  12  (1)  -- A  (SI,  32, - -  S  (I)  +1,...,S(N)) 

C  bii  A! — PRODUCT  OF  INDIVIDUAL  HEMS’  BOS  (I)  CB  ONI  DAI 

C  WHEN  ONE  IIEM'S  STOCK  LEVEL  INC  BIASES  BY  +1. 

C  BLUE  U,  J) — IN  CHASED  PCS  FOE  ITEM  I  WEEN  IIS  STOCK 
C  LEVEL  INCFEA5ES  EY  *1. 

C  B  K  A 1  1  ( J )  HAT  CN  DAY  J 

C  *.  PC  S —  WEIGHTED  VALUE  OF  SHAT  FOP  ONE  DAY 

C  A — SUE  CF  'WEIGHTED  POS <S)  OVEF  ALL  DC  AY  LAYS 

C  DEL  —  INCLEMENT  AL  INCREASE  IN  A  FOE  ONE  MCF.E 

C  IN I I  CE  ONE  ITEM. 

C  CElE AX--KAXIMUK  VALUE  OE  DEL 

C  CUE  (I ,  G ) — INDIVIDUAL  FKCEAEI1ITY  OF  S{I)  DEMANDS 
C  CV  IB  TEE  LEAD  TIME  FOE.  DAY  J. 

C  CUM  (I, J) —SUB  CF  EBCBAE1II11ES  CF  <  =  S  (I)  DEMANDS 
C  OVEF.  THE  LEAL  TIKE  FOE  DAY  J 

C  EX-- VALUE  CF  CUK(I,C)  FGF  CKE  DAY  A SC  CUE  ITEM 
C  PCUK--VALUL  OF  CUM  (1,1)  FCF.  CNF  DAY  AND  CNL  ITEM 
C 


**#*■+■#■*¥*■***#****■**** 


CCKMUN  EDDE ( 10  0)  ,SEX (100)  ,  EDCR  ST (100)  ,1 A MED  A  (100)  , 

*  V.  (ICC)  , ECCSI  (10  0)  ,SHCET  ,EAC  (10  0)  , 

♦COST, ECS (100) 

COMMON  DUiDAY , EUCGE1 ,TSHCE1 .TOC3IG  , 

*'C£1I  Cl. ,  TOTDKD  ,  E  FT  ,  K  EIGHT  ,  EEC  ,  S  Y5P0  S  ,  ND  A  Y  ,  NS 
COMMON  L L X  (ICC) 

CCKMOI.  ALAM  (20,  100) 

COMMON  Nil  Ei'.S  ,  NShUT  ,  SPIKE  { 10 C)  , ISP 
EFAI  r 

DO  U  El  I  PRECISION  DSDED,AS  (100)  ,  F.H  AT  B  ,  r  IDS  (S',  100) 

EC  0  E  L  £  PR  tCiSlOi;  B  b  AT  ,  V.PGS  r  DEL  ,DELM  A  X  ,  A  ,DE  I  ,  DELE  A  X 
DC  lib  LI  PBLCISIuN  L  OF.  (SO  ,  1 0  0 )  ,UJ  M  (2  0 , 1  2  C  ) ,  F  L  A  X  1  ( 1 C  0 ) 
EC  CELL  i?  BFC-iSICK  PX,PCUS 
COMMON  /STrilE/  MCAY(JC) 

iy£ni.D=1b3bE7.L^ 

ST  C  C  K  -  0  . 

C  G£’i  =  C  . 

a--:. 


c 


w  ®  *  *  v  *  Alii!  :M  ICE 

C 

DC  DU  1=  1,  HI  if  S 
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A2  (i)  =  :. 

CCS1=C0S1 +ECCST  (I) *LEX  (I) 

It  (BK1.EC.1)  GC  1C  52 

c 

C  *****BANDCXLY  DEIEKHIBE  KAIN1ENANCE  DAY  FCB  IT  EH  I 

C 

I.  =GGUDFS  (DSIED) 

JSP=IN1  (B*  (KDAY/2.-1)  +  .5)  +1 
M DAY  (I)  •=  JS  F 
SPIKE  (JSF) =1. 

KSr=JSP*ISi  (HDAY/2.) 

SPIES  (  P.  S  r )  -  1  • 

C  ALE  C  A  D  F  5  { -0 
SPIEL  <JE*)=C. 

SPIKE  ( K S F 1 =0. 

52  CO.iliKUE 

FCS  (I)  =0. 

SBa  (I)  =C. 

E  A  C  ( I )  =  C . 

50  CCNliKCE 
TSiiCBl-O. 

EEC=C. 

C 

C  *****  OPTION  =  2  KKEG  EVALUATING  EXIEI-.i.Alli  SImEIIAIL^ 

C  SEALES  MIX 

C 

IP  (OcT^O'S  k)  GC  10  105 
CCS1=0. 

DC  100  I—  1  r  Nil  Ef]5 
LEX  (I )  =  I AKEL  A  (1) 

COSI-COSI  +  0CQS1  (1)  *i-2X  (I) 

1U  CCtJllMJi 

105  COKllKUf 
1 1  AG- C . 

C 

c  *****  Bt GIN  AlGuailHK/LVAlUAllCl* 

C 

11C  DC  / - ~  J  =  1  t » D  •  •  Y 

lx  ( E I A u- 1 1.)  GO  10  265 
i.KEl-1 . 

C 

C  *****  FOB  EaCii  I  HE  Cii  DAY  J  COMFUTF  a  PCD 
C 

jO  250  1=  1  , lillEXS 

I  =  A I  A K (i,J) 

SAC  (I)  =d AC  ( I )  ♦  E D D l.  (I)  *EdCI.S1  (I ) /24. /ED ai 

i::  =  e}.f  (-i) 

t  c  l  h  =  I  X 
31  v_  i:.=  C  . 

II  (I  IX  (I)  .  Lx.  C)  GC  'I «.  7  ” 

I  =  I  E  X  ( I ) 

SlcC/:=0. 

L  C  leu  K  =  1  f  I. 

3ICCK=S1CC1>1. 

PX=FX*2/S!OCK 


1 


n  n  n  n  n  n  nnn  onnnno  onnn 


/ 
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PCJi1=PCUrl*EX 

BA  C  (I)  =  £  AC  (I)  -  (  1.-PC'J«)  /NDAY 
26  C  CONTINUE 

270  CONTINUE 

*****  CALCULATE  FBCDUC1  OF  INDIVIDUAL  ITEMS1 
FOB  CN  DAY  J 

Eh  AT=  FC  D  K*  E  H  AT 

PLUS  (I, J)-  (PCUM+  (PX*T) / (£TCCi>  1.)  )  /PCUM 
CUL  (1,  J)=FX 
cur,  (1,0)  =pcus 
ECS  JI)=PC’JK 

Si:GFT=  ( l.-PCUh)  *EDDE  (i) 

SEX  (I)  =5BX  (1)  +SHCET 
250  CCK12NUF 

BOAT!  (0)  =  Fri  AT 


*****  CALCULATE  WEIGHTED  PCS  (S)  CK  LAY  J 


WFCS=K  (J)  * E ii A 1 


*****  £ Ui'i  MIGLILL  PCS  FOB  All  DAYS 


A=litC£  +  A 

IF  (CE11GK.LQ.2)  GO  TO  281 
265  DEIKAX=C.C 

DO  280  1=1, NITERS 

*****  CALCULATE  A(£1,S2,...,S(I)  + 1  , .  .  .  ,  S  (r»)  ) 

EH  AT  2=EHAT  1  (J)  *PLUS  (I,  0) 

A2  (I) -A2  (I)  *  EH  AT  2**i  (J) 

IE  (J-LT.NDAY)  GO  TO  2«5 

*****  FILE  LMPhOVEKENT  OVEF  CONSENT  A(c) 


DZi=  (A2  (I)  -A)  / FCCSI  (1) 


*****  MM  OAX  iMPECVlhEIU 

IF  (LEL.LT.DELI1AX)  GO  TC  2c5 
L  E  E  5  T  =  j. 

EELKAX=LE1 

285  CONTINUE 

2  u  0  C  C  i>  T  I  .i  0  r 

2  Li  CONTINUE 


200  CONTINUE 

II-  (FLAG. 20. 1)  GC  TC  245 
CC  :  4,  1=1,!.  IT  ME 


c  *****  ADD  EES1  HEM  1C  SPARES  COST  IF  <  EUDGE1 

C 

1COS1-CCST  ♦  SCCST(IEEST) 

295  IF  (TCCS1.G1. EUDGZ2)  GL  10  530 
LCST=1C0S1 
A® A2  (IPES1) 

EC  300  1=1,N11£KS 
A2(I)=C. 

300  CGiiTiaVL 

C 

C  *****  hZL  EL3T  11 IK  10  SPARES  BIX 
C 

1EX  (IEEST) =  1EX (IEES1)  + 1. 

SUC3I=C. 

c 

C  *****  ELCC3PU1L  *E1GL1ED  PCS(S)  CN  DAY  0  FCF 
C  ME*  SPARES  fill  ARE  EACK08DEF  PROBABIII1IES 
.C  FOK  NEVLY  ADOBE  I1EE 
C 

DC  ECO  0*1, SEAY 
1=1EE£1 

Clili  (I,J)=CUR  (1/0)  *  ALAS  (J,J)/F1CAI  (LBX  (I)  ) 
RliAl  1  (0)=EHAT1  (0)  /CUM  (I,J) 

COE  (I,J)=CUK  (1,0)  +CUE  (I,  J) 

EAC  (1)  =EAC  (I)  -  (l.-CUK  (I,J)  )  /MAY 
PCS  (1;=C0B  (1,0) 

PH  All  (J>=EHAT1(0>  *CDM  (1,0) 

PiUS  (1,0)®  (CUB  (1,0)  +CUR  (1,0)*A1AK  (1,0) 
♦/FLOAT  (LBX  (I)  +1)  )  /COM  (I,  J) 

SHCR1=CUF  (1,0) *EDDE (I) ♦SBCR1 
500  CCNliHUI 

SEX  (I)=SEX  (I) -SIJGEI/WDAY 
ISIiCPI-ISbCFl-SHORl/NEAY 
CUBKl*NDAY-2+  (R£IGflT*2. ) 

C 

C  *****  CALCULATE  PCS  (S) 

C 

SYSPCS=A/COi-.»I 
E 1 A  3* 1 . 

C 

C  *****  FLAG  PBEVEK1S  FrCALCOLTlCS  OF  SE'a  EACKCRDSF 
C  FRCI ASillllLS. 

GC  TO  110 

530  IF  (SYSPCS.G1. .99)  GC  1C  54C 
C 

C  *****  IKCEEKEU1  EUDGII  CEILING 
C 

IF  (cDDGLl.Gl.lCClIG)  GC  1C  54C 
IUEGLl=EUnGET+1CC0000. 

LC  E31  1  =  1, MILKS 

EFr-=rIC*5AC  (J) 

—  31  C 1 1*  *  *  [!  E 

C 

C  *****  EiSPlAi  CCIPl'T 
C 


CALL  CUlrUI 
GC  1C  295 

545  CC  RC1— Kn^Y-*,+l.EI3H1*2 

£YSi?CS=A/CUE*l 

540  CCIiUN'JE 

DC  54c  I  =  1,\’12Efl£ 

EEO=EEC+PAC  (I) 

546  CCKLIKUE 
F  El  'JEN 
cl.D 
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c 
c 
c 
c 
c 

c 
c 
c 

C  THE  POLlOnliiG  V  AP1ABLE3  ARE  UNIQUE  T u  .10 LEI  R 

C 

C  CdG(K) — CD  KU-lAT  IV  E  WEIGHTED  PROBABILITY  CE  A 

C  H ACKGF.*>r L  1-CF  ITERATION  K 

C  ALPHA  —  .-tAXlliUK  ACCEPTABLE  PROBABILITY  OF  A 
C  BACKORDER  FOR  ANY  ITEK 

C  ?.’*EPS —  NUMBER.  CE  ITERATIONS  DESIRED 
C  MAX  —  MAXIMUM  IT  Eh  COST 

C  C'JfHT — CUMULATIVE  WEIGHTS  FCF.  L ACKCF DELS  ON  ALL  DAYS 
C  COJSPAh  —  ID DICaLiiS  SUBROUTINE  IS  BEINu  USED  TO 
C  EVALUATE  EXTERNALLY  GENERATED  SEAFLS  MIX 

C  EXBO  (L  ,  *.) — EXPECTED  WEIGHTED  EACH  CEDE  RS  CP  ITEM  I  AT 
C  IT  IF  ALIGN  K 

C  IS1CCK  il  ,K)  --S1CCK  LEVEL  FCR  ITEK  I  AT  ITERATION 
C  CCSTM  ((■.) — SPA? IS  KiX  C0S1  FOE  ITERATION  K 

C 


SUE F0U1  IKE  TO  IKPI  EMENT  ALGORITHM  ALL  EVALUATE  SPARES 
MIXES  i-CL  HOLE!  p 

EY  KATHLEEN  CCNLEY,  DECEMBER  1SB1 

SUEEOUIINE  EE 


COSMOS  BDDR  (  1 G 0) #FDCRS1 (100)  , 

*  W  (ICC)  , HCCST  (IOC)  ,EG  <30,15)  ,CEC  (15)  , 

*CCST  (2) 

CCXJOl  IX  ,  SUED  AY , ALPHA , N  REPS , EJDSET  ,  K  AX  ,  CUMUl , 
* CPT 1CK,  T OT  DM C,  CCKP  Ah,  WEIGHT,  DDF.  ,NDAY,iJS 
COMMON  LEX  (1  CC)  ,CUKBC  (20)  ,EXEC  <3C ,2C) 

CC  K  HO  S  AL.AM  (100,100)  ,L  STOCK  (30,20)  ,  COS  IK  (20) 

C CUM ON  SI  TEDS, NSHDT, SPIKE  (ICO)  ,ISP 
REAL  r-X,  PCUM,  THETA,  ThETAK 
DIMES SICS  P  (ICO) 


C  *****  ELSE  MAX  COST 

C 

IF  (CPI  ION. NZ. 3)  GC  TO  215 
MAX=C. 

DC  Jou  I=1,i;iTEKS 

Ii:  (HCCST  (I)  .1E.KAX)  GC  TC  25C 
i.  AX  =  HCCST  (I) 

2  5 C  couth;  UL 

WRITE  (50,305)  (ALAR  (I,J)  #J=1,UD;,i) 
305  FCF MAT  (15  (/7F10.5) ) 

300  cc  m  i  ;;u  z 
c 

C  *****  FILL  sor:  H  t. ZIGIiT s 

c 

c  u t .  —  j  m 

DO  J  ~'  j  J  *  1  ,  i*  */  rt  i 

cut’.n-cuHv.T+v:  (J) 


U  (J  U  U  UUU  OL>U 
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353  CONTINUE 
315  CCMIulli 

TiiEl  A  =  (1  .-ALPHA)  *CUKb1/EAX 
EC  355  M  M  =  1 ,  N  t  E  F  £ 

CUMBC (ME) =0. 

CCSTM (MM) =0. 

C BO  (ME)  =0 
355  CCR1IXUE 

DC  5CG  1-1  ,  N 1 1 E  N  S 
CUXSUM=G. 

S1GCK=G. 

N=KPIFS 

It  (COME  AS. Eg. 1)  N= 1 
EC  4  50  M=1,:v 

*****  COMFU1E  MU  LI  IFL1EE  FOB  1IEPATION  ti  AND  CRITICAL 
V  A 1 U  £,  EOS  11  EK  I. 

1  HL1  AM  =  Iiit;i  A/2.  **  (K-1.  ) 
SHS=CUMft’T-lHilAM*BCOST  (1) 

It  (51CCK.G1.0)  GC  10  3 b0 
CUhSUM=0. 

*****  S’JM  FBOE AE111IES  OF  ZEBC  LEAD  TIKE  DEMANDS 


37  0 
360 


35  0 


370 
_>7  7 


EC  370  J=  1 , N  D  AY 
1  =  AL AK  (I,G) 
t  (J)  =  EXP  (-T) 

ECUH=E  (J) 

CUSSUM-CUMSUM  +  ii  (J)  *FCUM 
CC Nil  NOE 
CCMINUE 

IF  (CCHEAB. EQ. 0.)  GO  IO  3S0 
EUS=C. 

SICCK=LBX (I) -1. 

GC  1C  355 
CONTINUE 

IF  (CENSUS. G£. BBS)  GC  1C  4CC 
C  H  M  5  U  M  =  0  . 

S1CCK=STCCF+1. 

EC  375  J=1,NDAY 
1  =  ALAK  (I  ,  J) 

F  (J)=EXE>  (-T) 

EC UM=  E  (J) 

IF (STOCK. EQ.O.)  GC  TC  377 
I  SI CC  K= INI  (STOCK) 

DO  376  K=1,IS10Cr 

P  ( J)  =F  (J)  *1/  FLOAT  (K) 
FCUK=F  (J)  +  FCLV. 

CCfllKUE 
C  C  .i  1 1.  N  U  L 

v.  — .J  v  IS  OF  4-  •.  i  1  iiCh  Ci>  no 

CUMSUK=CUN.SUti*K  (J)  *FCUM 
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I 

L 


375  CONTINUE 

C 

C  *****  FZfZAT  UNTIL  STOCK  CAUSES  FUNCTION  TO  ECXC2E0 
C  CK11ICAL  VALUE. 

C 

IF  (CCKPAB.EC.  0. )  GO  TC  380 
400  CONTINUE 

EXEC  <1  f  S)=0. 

1STCCK  =  S 1 OC K 

COSTS  (M)=CCSTK  (K)  ♦  LSTOCK  *ECOST  (1) 

C 

c  *****  EVALUATE  LXtxCTED  SLIGHT EU  EACKC&LFES  FCH 
C  EACH  IT  ES 

C 

UO  4  60  J=1#NDAY 
3UK=C 
ST  K=ETOOK 
EK=P  |J) 

T  =  ALAS  (I, J) 
ma**  STK=Sli;*1. 

ib=ih*l/S1K 

S'JM=SUN  ♦  (S1K-STOCK)  *Pii 
IF  (PK. IT. IE-5)  GC  TO  465 
GO  TO  4b2 

465  CONTINUE 

EXEC  (I, is)  =EXEO  (I,h)  ♦  *  (J)  *SUK 
46  C  CONTINUE 

C 

C  *****  EVALUATE  CUMULATIVE  Ej.  ST  EC] 

C 

CEO (£) =CEO  (H) +  EXEC  lI,a)/CU KH1 
450  CONTINUE 

50  C  CONTINUE 
CALL  CUi 
KETUriN 
EN  E 
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SUoZiCUTIiiL  10  IKPLjiKSNI  AlGCtllhK  AND  EVALUATE  SPALLS 
MIXES  FOR  MODEL  D 

EY  K  AlULtEN  C C i;  1 E  Y  ,  DECEhEEE  19a1 

SUEROUTIRE  DD 


C  THE  FOLLOWING  VAFIAPLES  A  EE  DEFINED  E  OK  NEITHER  ftOEnL 
C  A  ft  Cn  MODEL  £: 

C 

C  ,.EAF  (I,  J)  — WIIGliT  II  A  EACKCEDEP  DAY  UHLS  EXPECUD 
C  WAITING  TIMS  E  Op  A  E ACKCRDZE  o.i  L  AY  J  FOE  I1IL 

C  I. 

C  CRT  (1) — CUMULATIVE  REAR  {1,0}  ICR  ALL  I«YS  cOk  II  EE  T. 

C 

C*********** ************** ***************************M******** 

c 

CCKttON  BDDB  (100}  ,RDCEST  (ICO)  , 

*  ri  (100)  , ECGS1  (100)  ,EO  (3  0,  15)  ,CBO  (15), 

*CC£i  (3) 

CCKhOK  PX, CO  ELAY , ALPhA, KEEPS, BUDGE!, MAXI , MAX , 

♦OPTION  ,  IClDEiD  ,CGr;fcrtR,KiilGiiI  ,  DDfc  ,  ND  AY  ,  NS 
COMMON  LEE  (  100)  , CUMBC  (20)  ,  EXEC  (30,20) 

CO ft MOD  ALAK  (ICC, ICO)  ,LS1CCK  (30,20)  ,CCSlft  (20) 

COMMON  NITEMS,  NSHUT,  SPIKE  { 100)  ,  ISP,  LEAF.  (30,100)  ,CWT  (30) 
REAL  P X  , PC  Uft , ‘i HE! A  ,  TilEl  AM 
DIMENSION  E(1GC) 


IF  (OPTION. N£. 3)  GC  1C  315 
C 

C  *****  FIND  'SUM  OF  WEIGHTS 

C 

CLKK1=0. 

DO  3  50  O  —  1 ,  N  £  AY 

CUMKl=CDUkl*N  (J) 

350  CONTINUE 

L*  O  JLo  I =  1  ,  I  L  L  w 
C 

C  *****  COMPUTE  InEAD  (1,0) 

c 

C  Vi  1  (1)=C. 

DO  3  10  J=  l,K£AY 

t:  BAR  (1 , 0 )  =  K  (J)*E£CRS1  (I)  /24.  /  AL  AM  (I,  J)  /COMST 
CWT  (I)  =Ca1  (I)  ♦  is  CAR  (I,J) 

310  CONTINUE 

2  RITE  (50,  100)  C*1  (I)  ,1,  (NEAR  (1,0)  ,  J  =  1  ,NDAY) 

ICO  FCEHAI  (E1C.  C, 14,15  (/7F10.5)  ) 

300  CCNIaNUL 

c 

Y  I;  I Y  A  =  1  0  0  0  o  U  0  0  0  * 

DC  11C  1=1, LITERS 


no 


C  *****  FIND  MAXIMUM  TtiETA  FOB  AIL  HEMS 
C 

1  d  El Al  =  (1-A1FBA) *Cfil  (D/ECCSI  (I) 
IF  (THLTA.GI.TbcTAI)  1 HET A=THETAI 


110 

CONTINUE 

WRITE  (60,  105)  THETA 

105 

FORMAT  (• IK£TA=  •  ,F10 

315 

CONTINUE 

DC  355  MM= 1 ,  NBEES 

CU MBC  (MM) =0. 

CCSTM  (MM;  =0. 

CEO  (MM)  =0 

355  CONTINUE 


EC  500  1=1, NIT EK5 
CUM5UM=C. 

STOCK=Q. 

N=NR£?3 

IF  (CCMFAF.EQ.  J)  N=  1 
DC  45 C  K=1,K 

THET AM=1HETA/1. **  (M-  1.) 

BhS=CKl (II -T  H  El  A  M*  EC  CS  T  (I) 

IF  (STCCI'.GT.O)  GO  ‘10  3d0 
CUKSUM=0. 

EC  37 C  J=1,NDAE 
C 

C  *****  SUM  FS03AI31LIIIES  Cr  ZEBC  DEMANDS  OVER  LEAD  TIME 

C 


non  non  nnon 
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ICUK=E  (J)  +  PC  Us, 

370  CONTINUE 

377  CONTINUE 

cahsufl=cunsa«*k’3is  (i,  J)  *pcus 

375  CONTINUE 

If  (CCRPAfi.EQ.  C.)  GO  TO  380 

*****  EEPEAT  UNTIL  STOCK  CAUSES  FUNCTION  TO  EXCEED  CEITICA1 
VALUE 


4CC  CONTINUE 

ESTCCN  (I, SI)  = STOCK 

CCSTE  (K)=CCSTM  ( M )  +1S1CCK  (I,«)  *ECOST  (I) 

I X  £  C  (i, h) =0. 

*****  EVALUATE  EXPECTED  WEIGHTED  BACKCuDEtf  DATS  FCE  EACH  ITEM 

DC  460  J=1,NBAX 
S  U  2=  0 
ST K= STOCK 
1  =  ALAK  (I, J) 

P?=P(J) 

ST  F-ST  K+  1 . 

Ff,  =  Fr*T/STK 

SUK=SUM+  (STF-STCCK) *  E  h 
IF (PL. LT.  11-8)  SC  TO  4o5 
GO  TC  462 
CONTINUE 

EXEC  (I  rM)=EXEC  (I,K)  ♦WEAL  (I,  J)  *SUa/NDAY 
CONTINUE 

*****  EVALUATE  CUMULATIVE  EXPECTED  WEIGHTED  B ACKOLDEt  DAYS 

CUKEC  (K) =CUMEO (K)  +EXBO  (1,55) 

C DC  (?)  =CU 53 EC  (E) 

450  CONTINUE 

53 C  CONTINUE 
CALI  OUT 
R  £ 1  U 1.  K 
ENL 


4bz 

46  5 
46  0 
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c 

C  SIHU1ATI0N  fOE  EASING  STUDY  FOB  SHUTTLE  SPARES 
C  BY  KATHLEEN  CCNLEY,  DECEREER  1981 

C 

C  HOSE:  SHE  RAIN  PBOGBAB  IS  ALHOST  ENTIRELY  TAKEN  FBOH  THE 
C  BEFEBENCE  GIVEN  BELOW .  IT  IS  INCLUDED  FOB  PUBFOSES 

C  OF  VARIABLE  DEFINITION  ONLY. 

C 

C  BEFEBENCE:  SIROLATION  OF  A  TKO-ECHEION  INVENTORY  SISTER 

C  BY  JIR  COGL1 ANO ,  NOVERBEB  1980 

C 

C*********************** ************************************** 

c 

C  LIST  OF  GLOBAL  VABIABLES 
C 

C  / PAEAR/  SISTER  FAiARETEBS 

C  LHAZ  NOHBEB  CF  BASES 

C  IR AZ  NUHBEE  Of  ITER  TYFES 

C  N UNITS  (L , I)  NO.  Cf  ONUS  CF  TYFE-I  AT  EASE-L 

C  NSEABE  (L, I)  INITIAL  NO.  OF  SPARE  PARTS  OF  TYPE-I  AT  L0CATI08-L 
C  AFAIL  (L,I, J) AVG.  TIHE-TO-F AILOBE  FOB  TYFE-I  PARTS  INSTALLED  IN 
C  TYFE-J  ONUS  AT  BASE-L 

C  AREP  (1,1)  AVG.  REPAIR  TIRE  FOR  TYFE-I  PARTS  REPAIRED  AT  LOCATION-L 

C  PREP  (1,1)  PRCB.  OF  A  BASE  BEPAIB  FOB  TYFE-I  PARTS  FAILING  AT  BASE-L 

C  ABD(L)  AVG.  TRANSIT  TIRE  FEOR  BASE-L  TC  THE  DEPOT 

C  ADE(L)  AVG.  TRANSIT  Tint  FP.Ofl  lilt.  utFCT  TO  BASE-1 

C  AORD(L)  AVG.  ERCCESSING  TIRE  FOE  AN  CBEEfi  FBOR  BASE-L 

C  1DEPOT  1NDEZ  FOB  THE  DEPOT  (ALWAYS  EQUALS  LRAZ+1) 

C 

C  /STATE/  STATE  VABIABLES 

C  NB  (L , I)  NO.  OF  TIFE-I  PABTS  IN  BEPAIB  AT  LOCATION-L 

C  NS  (L, I)  NO.  OF  TYPE-I  PARTS  IN  SPARE  STOCK  AT  LOCATION-L 

C  NBD  (L ,1)  NO.  OF  TYPE-I  PARIS  IN  TRANSIT  FBOH  BASE-L  TO  THE  DEPOT 

C  NDB  (1,1)  NO.  OF  TYPE-I  PARTS  IN  TRANSIT  FBOH  THE  DEPOT  ^0  EASE-L 

C  NU(L,I,J)  NO.  CF  TYPE-I  PARTS  INSTALLED  IN  TYPE-J  UNITS  AT  EASE-L 

C  NG  (L, I)  NO.  Of  TYPE-1  UNITS  AT  BASE-L  GROUNDED  FOR  LACK  OF  FARTS 

C  NGZX  NO.  GROUNDED  FOB  LACK  OF  FARTS  OVER  ALL  UNITS  AND  BASES 

C 

C  /STATS/  STATISTICS:  CUHULATIVE  SUttS  Of  I  BE  STATE  VABIABLES 

C 

C  /CLOCKS/  SIRUI AT1CN  TIRING  VABIABLES 

C  CLOCK  CURRENT  TIRE 

C  CST ART  LENGTH  CF  HABR-UP  PERIOD 

C  CSTOF  TIHE  TO  STOP  SIHULATIUN 

C  CINTES  LENGTH  CF  TIRE  BE1EEEN  STATISTICS  OBSERVATIONS 

C  HOES  NONEEB  GE  OBSERVATIONS  TO  TAKE 

C 

C  /OPTION/  OPTIONAL  FEATURES 

C 

C  IV IB  1  FOB  VARIANCE  ESXIRAIION 

C  I EXP  1  FOR  EXPONENTIAL  FAILURE  TIRES  (FASTER  EXECUTION) 

C  HR  ACE  1  FOR  TRACE  OF  EVENTS 

C 

C  /POLICY/  P01ICY  SELECTION  VABIABLES 

C  IBOBD  '  BASE  ORDERING  POLICY 


U<JOU<JU<JUOUU  uo  u«~ 
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C  1BINST  BASE  INSTALLATION  POLICY 

C  ID5BIP  DEPOT  SHIPPING  POLICY 

C 

C  /OBCEBS/  DEPOT  BACK-OBDEB  LIST  ( ABBAZ  OF  FOBNABD-LINRED  LISTS) 

C  /EVENTS/  FUTOBE  EVENTS  LIST  (HEAP-S0B1 EE) 

C  /SEEDS/  BAHDOB  NUNBEB  SEEDS  (DOUBLE  PBECISION  FOB  IBSL) 

C  /CUBBEV/  CDBBENT  EVENT  CODE  AND  ATTBIEUTES 

C  /NB/  NEED-BELOCTAHCE  V ABIAELES 

C  /FILES/  INPUT  AND  OUTPUT  FILE  UNIT  NUHEEBS 

C 

C*m»**M*Mm*nM*M»*MM***M**M**M*MM*m*M*M»m 

c 

COMMON  /PAfiAM/  LMAX,  WAX,  NUNITS(5,5),  NSPABE(6,5), 

5  AVAIL  (5,5,5)  ,  ABEP(6,5),  PBEP(5,5), 

6  AED (5)  ,  ADE  (5)  ,  AOBD  (5) ,  LCEPOI 

COMMON  /STATE/  NB(6,5),  NS(6,5),  NBD(5,5),  NDB  (5,5)  ,  NU(5,5,5), 

S  NG  (5,5) ,  NGXX 

COMMON  /STATS/  SB  (6, 5),  SS(6,5),  SBD(5,5)#  SDE(5,5),  SO  (5,5,5), 

S  SG  (5,5) ,  SGXX,  SSGXX, 

5  SSB  (6,5) ,  SSS (6,5)  ,  SSBD(5,5),  SSDB  (5,5) , 

£  SSO  (5,5,5) ,  SSG  (5 , 5) 

COMMON  /CLOCKS/  CLOCK,  CSTABT,  CSTOP,  CINTEB,  NOBS 
COMMON  /OPTION/  IVAfi,  I2XP,  1TB ACE 
COMMON  /PCLICZ/  1E0BD,  IBXNST,  IDSBIP 

COMMON  /OHDEBS/  ITPBH  (5,2) ,ITPBT  (5,2) ,10C (100) ,NXI10C(100) , NXTLF 
COMMON  /EVENTS/  NXTEV(lfinO),  T^MIOeO)  .  INFO  (1000,4)  ,  SXTEVF, 

C  FTIHE,  FBATE 

COMMON  /COBBEV/  NODE,  K1,  K2,  K3 
COMMON  /NB/  BNMAX  (5) ,  LENBAX  (5)  ,  EN(5,5), 

6  B1AELE (5, 5, 2 1 ,2 1) ,  DT ABLE (21) ,  NIHAX,  NSHAX, 

6  ICB  (5,8,4,4,4,4,4) ,IDHAX,  IEMAX 

CCHMON  /FILES/  JIN,  JOUI 
COMMON  /LATS/  S1,S2 


***** 

TBE  MAIN  PBOGBAB  CONTEOLS  THE  EXSCDT1CN  OF  TBE  SIMULATION. 

THEBE  ABE  THBEE  PHASES: 

(1)  BEAD  THE  INPUT  PABAMETEBS  AND  INITIALIZE  THE  ST  ATE -VABIABLES. 

(2)  PBOCESS  THE  EVENTS  AS  THEY  OCCUR.  ADVANCE  THE  CLOCK,  COLLECT 
STATISTICS,  AND  CALL  A  ECUT2NE  TC  CABBY  0D1  TBE  DETAILS. 

(3)  PBINT  THE  SUMHABZ  STATISTICS. 

***** 

JIN  =  1 
JOOT  *  2 


CALL  INPUT 
CALL  IN1T 
CALL  ST  AT  (1) 

CALL  EVENT  (T,K0CE,K1,K2,K3) 
CLOCK  *  T 
C 
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IF  (K0DE.EQ.1)  CALL  AEBIVE  (K1,K2,K3) 

XI  (K0DE.EC.6)  CALL  FALL  (K1,K2,K3) 

C 

C  *****  LAUNCH  CBITICAL  EVENT 
C 

If  (XODE.EQ. 15)  CALL  LCHIT  (K1,K2) 

If  (NODE. Eg. 15)  CALL  CBDEE  (K1,K2,K3) 

If  (NODE. EC. 16)  CALL  BEFAIB  (K1,K2) 

Xf  (NODE. EC. IS)  CALL  STAY  (2) 

If  (NODE. N£. 19.CB.CL0CK.L1.CS10P-C1NTEB/2)  GOTO  10 
C 

CALL  ST AT  (3) 

C 

STOP 

BBS 

C 

c  *****  COBPOTEB  CODE  FOB  ALL  SOI  THE  LAUNCH  CHITICAL 
C  SDBEOUTINE  IS  OBITTED--HOHEVEB,  TO  BUN  THE  PfiOGBAB , 

C  ALL  THE  EVENTS  BUST  COBBESfCNC  TO  SUEEOUTXNES 
C 

c 

c********************** *************************** ***•••••**** 
c 

C  LAUNCH  CRITICAL  AND  NON ST AT ION ABY  DEBAND  SUBBOUTINE: 

C 

C«**M***MM»*M**MM«MMM*M*M»^**M»MM*»**»MM**M* 

c 

SUBBOUTINE  LCBXT  (L,K) 

C 

COBHON  /CLOCKS/  CLOCK,  CSTABT,  CSIOf,  CINTEB,  NOBS 

COBBON  /EVENTS/  NXTEV  (1000)  ,  TIHE(ICOO),  INFO  (1000,(1)  ,  NXTEVF, 


C  FTIBE,  f BATE 

COBBON  /PABAB/  LBAX,  IBAX,  NUNITS(5,5),  NSPABE(6,5), 
t  Af AIL  (5,5, 5) ,  ABEP  (6,5)  ,  PBEP(5,5), 

6  ABO  (5),  ADE  (5) ,  AOBD  (5) ,  LEEPOT 

CCHBON  /STATE/  NB(6,5).  NS(6,5),  NBD(5,5),  NDB(5,5),  N0(5,5,5) 
G  NG  (5,5) ,  NGXX 

•  COBBON  /POLICY/  IBOBD,  IBINST,  IDSHIP 
COBBON  /COBBEV/  KCDE.K1 ,K2,K3 
C 

C  *****  R  is  1  BHEN  LAUNCH  CBIIICAL,  BESET  TO  0  UPON  LAUNCH, 


C  AND  IS  2  BHEN  LAUNCH  CBIIICAL  IS  APPBOACHING 
C 

c  *****  EETEBSINE  BASE  TO  BASE  TBAVEL  IIHE 
C 

LL-L+2 

BS1=ADB  (LL)  •  2 
If  (ADB(L).GT.BST)  BST=ADB  (L) 

C 

IFJK.EQ.  1)  GO  TO  100 

c 

C  *****  CUEBY  USEE  FOB  NONSTATION ABITY  IN  CEflAND  SATE 
C 

IF  (TFLAG.GT.O)  GO  1011 
HBITE  (7,9) 
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S  FCE.5A1  {‘IIKES?  F  3 .  C '  ) 

EEAt(o,10)  TIKES 
1C  FCFMA1  { t  3 . C) 

11  I F 1 AG=  1 . 

C 

IF  (ft. tg.2)  GO  IG  30 

c 

C  *****  If  iJUJNCfj  IS  TODAY  (K  =  0)  ,  DOUBLE 
C  I  Hi  EEMANC  EA1F  FCF  TODAY  ONLY 
C 

FT 1K£=FT1XE -TIKES 

IF  (I  Ti.'iE.  IT.  C)  1-TIM=C. 

C 

C  *****  SCHEDULE  ANCluZi  LAUNCH  CEI1ICA1 
C 

T=  5C/K  UNITS  (1  , I K  AX )  -  EST 
LI  Shir  =  L 

CALL  SCliHO  *1,12,1,1,0) 

C 

c  *****  sciiLLlJIx.  ANC1HEF  LAUNCH  CFI1ICAL  AFFEC  ACliING 

C 

12  =  1- ADE  (II) 

IF  (T 2-  GL.  C)  GC  TC  25 
SEITZ  (GGU'I ,  2  C ) 

20  FCE'SVI  ('LAUNCHES  ICO  CLOSE  FOE  LA1ILA1  r.ESUrFlY') 
SICE 

25  CONTINUE 

C All  SCHLD  (12,12,1,2,:) 

C 

C  *****  It  LAUNCH  CblllC  AL  A F F EC  AC HI N3 , PCSSI E L Y  SCUEEUL 
C  ShlFXENT  TC  THE  EASE 

C 

30  DC  50  1=1,  IK  AX 

I S  ri  I  ?  =  I 

ii  3=1 

Call  L5NIF  (ISHIE  ,ISHIF,  JCEL) 

CALI  CACT  (IShlF , ISnlr , JOED) 

5C  CONTI  Mis 

Y.  3=0 

GC  TO  y  'j  k, 

C 

c  *****  t:  r r iiSLNILY  IAUNCH  CFITICAL,  SChiEUIE 
C  A  LAUNCH  ( K  =  0 )  IN  EST  TIKE  UNITS 
C 


